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Abstract  of  the  Dissertation 

Performance  Feedback  &  Control  of  Solar  Concentrators  using 
Wave  Front  Sensing  Techniques 
By 

Joseph  N.  Beasley 

Claremont  Graduate  University:  2007 


The  major  requirement  for  using  concentrating  Solar 
Thermal  devices  is  the  proper  placement  of  the  focal  spot 
on  the  absorber  to  provide  heating  of  the  working  fluid  to 
produce  thrust  or  to  generate  electricity.  Without  proper 
placement  of  the  focal  spot,  either  in  a  spacecraft  or 
ground  unit,  solar  energy  is  not  transferred  to  the 
propellant  or  working  fluid  or  at  best  case,  a 
significantly  smaller  proportion  of  the  incident  energy  is 
transferred  to  the  fluid. 

A  control  system  is  needed  for  ground  units  to  track 
the  Sun  across  the  sky  or  to  track  the  Sun  as  the 
spacecraft  maneuvers.  The  control  system  would  use  a  fine 
focus  controller  after  acquiring  a  general  location  for  the 
Sun.  Optimal  focus  was  defined  for  a  space-based  system  as 
0.1  degree  and  0.1  inch  in  x,  y,  and  z. 
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The  research  problem  is  to:  (1)  develop  a  sensor  that 
provides  information  about  the  location  of  the  focus  beam 
from  the  concentrator  that  also  can  tolerate  the  extreme 
temperatures  involved  in  providing  solar  heating,  and  (2) 
develop  the  mathematics  of  a  control  system  that  will  focus 
the  concentrator  in  an  optimal  fashion  based  on  the  sensor 
measurements.  The  sensor  developed  in  this  dissertation 
solves  the  temperature  issues  by  using  the  absorber  as  a 
wave  front  Sensor  and  imaging  that  device  with  a  Charge 
Coupled  Device  (CCD)  to  determine  the  location  of  the  focal 
beam . 

A  hardware  model  of  the  absorber  and  concentrator  was 
used  to  generate  images  for  analysis.  Images  generated 
with  the  hardware  models  are  analyzed  by  first  searching 
the  images  for  intensity  areas  on  the  sensor,  and  then 
profiles  of  the  intensity  areas  are  integrated  to  find 
centroids  of  the  intensity  areas.  Once  the  centroid  is 
found  for  the  intensity  areas,  the  controller  can  generate 
commands  to  move  the  focal  beam  to  proper  focus. 

Results  from  this  research  indicate  that  the  modified 
wave  front  Sensor  method  worked  and  provides  control 
information  to  properly  control  the  concentrators  to 
provide  optimal  heat  transfer  to  the  heating  fluid. 
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INTRODUCTION 


For  direct  thrust  solar  thermal  propulsion  or 
terrestrial  solar  heating,  proper  positioning  of  the  solar 
focal  spot  is  paramount.  Along  with  positioning 
information,  a  sensor  system  must  tolerate  the  intense  heat 
generated  by  a  solar  concentrator.  The  sensor  system 
developed  in  this  dissertation  satisfies  both  of  the 
previous  requirements.  The  system  developed  in  this 
dissertation  utilizes  the  thruster  -  absorber  as  the  focal 
spot  sensor  with  a  Charge  Coupled  Device  (CCD)  camera  to 
image  the  absorber  -  sensor.  Two  image  processing  methods 
are  applied  to  the  CCD  images  to  generate  the  focal  spot 
location.  Once  the  focal  spot  is  located  on  an  image 
positioning  information  and  commands  are  generated  to 
reposition  the  concentrators.  Both  methods  produce  the 
focal  spot  information  needed  to  generate  positioning 
information  for  the  concentrator.  This  dissertation  will 
follow  the  experiments,  research,  and  analysis  utilized  for 
the  development  of  the  new  solar  sensor  system.  It  also 
introduces  the  term  "blobs"  to  refer  to  intensity  areas  of 
light  reflected  onto  the  tubes  of  the  sensor  from  the 
concentrator.  These  blobs  are  images  of  the  light  in  the 
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concentrator  and  also  are  related  to  the  intensity  or  power 
being  reflected  to  the  sensor. 

Chapter  1  Solar  Thermal  Propulsion  Back  Ground 

The  solar  concentrator  field  may  be  divided  into  two 
primary  areas  of  interest:  terrestrial-based  and  space- 
based.  Both  types  of  concentrator  can  be  either  imaging  or 
non-imaging  and  they  can  be  rigid  or  inflatable. 

Some  of  the  applications  of  solar  concentrators 
include:  terrestrial  power,  space  power,  material 
processing,  hazardous  material  destruction,  building 
illumination,  space  plant  illumination,  and  other 
applications . 

Solar  Thermal  Power  as  propulsion  has  its  roots  in  the 
work  of  Krafft  Ericke's  Solar  Powered  Spaceship  done  in 
1956.  Ericke's  spaceship  was  the  first  design  that 
utilized  a  double  solar  concentrator  connected  together  by 
a  hydrogen  powered  thruster.  Further  development  of  the 
solar  powered  spaceship  concept  was  sporadic  during  the 
1950's  and  1960's  until  a  single  concentrator  thruster  was 
developed  in  1963.  After  1963,  the  next  development  of  the 
technology  occurred  in  1979  with  the  development  of  a  twin 
off  axis  paraboloid  concentrator  concept,  by  Rocketdyne. 
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The  solar  thermal  field  was  then  relatively  dormant  from 
1979  until  the  Rocketdyne  thruster  tests  of  1984-1991  with 
a  new  Rhenium  coil  cylindrical  cavity. 

The  solar  thermal  propulsion  activity  really  took  off  in 
the  1990's  with  the  development  of  castable  thin  films  for 
casting  inflatable  concentrators.  The  1990' s  also  saw  the 
development  of  rigidized  structures  for  constructing  the 
support  structures  for  inflatable  solar  concentrators.  The 
year  1997  saw  the  successful  deployment  of  a  2  X  3  meter 
inflatable  concentrator  that  signaled  that  the  technology 
was  ready  for  up  scaling.  The  development  activity  from 
1997  onward  has  included  the  manufacture  of  a  mandrel  for 
casting  a  4  X  6  meter  inflatable  concentrator. 

Solar  Thermal  Propulsion  (STP)  is  a  promising  concept 
for  use  in  an  Orbital  Transport  Mission,  a  deep  space 
mission,  or  other  upper  stage  space  missions  requiring  high 
Specific  Impulse  engines.  Specific  Impulse  or  Isp  is  a 
rocket  propulsion  term  that  measures  the  thrust  per  unit  of 
propellant  expended.  Thus  STP  systems  have  a  much  higher 
Isp  than  chemical  systems  of  the  same  size.  Since  the  energy 
for  propulsion  is  available  to  the  spacecraft  in  orbit,  in 
the  form  of  sunlight,  an  STP  system  would  not  have  to  carry 
both  propellant  and  oxidizer  to  produce  energy  for  the 
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upper  stage  unit.  Only  an  atomic  or  molecular  propellant 
such  as  hydrogen  would  be  needed  on  the  upper  stage. 

The  STP  system  heats  up  hydrogen  in  the 
thruster/absorber  and  then  expands  the  hydrogen  through  a 
nozzle  to  produce  the  thrust  required  by  the  vehicle.  Thus 
the  thrust-to-weight (of  propellant)  efficiency  of  the  STP 
system  is  much  better  (more  than  two  times  better)  than  the 
equivalent  chemical  upper  stage  as  less  non-payload  mass  is 
needed  in  the  STP  system  for  the  same  amount  of  thrust  to 
be  developed.1 

The  major  difference  between  a  chemical-type  thruster 
or  spacecraft  and  an  STP  thruster  or  spacecraft  is  that  the 
STP  spacecraft  does  not  carry  the  added  weight  of  the 
oxidizer  as  the  chemical  spacecraft.  Thus,  the  STP  has  a 
weight  advantage  over  the  chemical  spacecraft  in  that  a 
lighter  molecular  mass  propellant  is  used. 

The  STP  system  also  has  the  added  task  of  controlling 
the  solar  concentrators.  Not  only  does  the  spacecraft  need 
to  properly  point  the  concentrators  towards  the  Sun  while 
the  craft  maneuvers  in  space,  it  needs  to  protect  them  from 
the  exhaust  of  the  spacecraft.  These  two  requirements 
limit  the  control  region  of  the  concentrators,  thus 
possibly  limiting  the  available  power  or  limiting 
maneuverability  of  the  spacecraft  itself.  A  final  control 
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requirement  is  that  the  focal  cone  should  always  be 
positioned  on  the  aperture  closest  to  the  nozzle  for 
effective  heat  transfer  to  the  propellant.11  The  focal  cone 
is  defined  as  the  envelope  of  concentrated  light  that 
approximates  the  hyperbola  of  rotation,  whose  surface 
intersects  the  concentrator  perimeter  and  the  boundary  of 
the  region  of  maximum  intensity.  The  focal  cone  changes 
shape  and  position  when  the  concentrator  is  shifted  from 
the  optimal  or  focused  position  see  Figure  1. 


Figure  1  Hyperboloid  Construction 

The  current  AFRL  concept  for  the  concentrators  for  a 

solar  thruster  is  two  off  axis  paraboloid  concentrators 
connected  to  a  central  thruster.  See  Figure  2. 
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Figure  2  STP  SPACECRAFT  CONFIGURATION 


Each  concentrator  is  connected  to  the  central  thruster  with 
a  hexapod  or  Stewart  platform  offering  6  degree  of  freedom 
(DOF)  control.  The  6  DOF  needed  are  yaw,  pitch,  roll,  x, 
y,  and  z.  The  axes  and  angles  are  depicted  on  the 
schematic  shown  in  Figure  3. 


Figure  3  6  DOF  Axes  For  Spacecraft 
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Optical  analysis  for  a  flight  demo  concentrator  system 
indicated  that  the  required  arbitrary  angular  pointing 
accuracy  is  0.1  deg  and  the  translation  accuracy  should  be 
0.1  inch.  Coarse  alignment  is  obtained  using  an  on  axis 
detector,  such  as  a  Sun  tracker  or  photo  diode,  which 
brings  the  spacecraft  and  concentrators  into  a  general 
alignment  with  the  Sun.  Coarse  alignment,  within  a  few 
degrees,  brings  the  focal  cone  into  rough  alignment  with 
the  desired  aperture  on  the  thruster.111 

A  spacecraft's  coarse  alignment  brings  the  focal  cone 
into  rough  alignment  with  the  desired  aperture  on  the 
thruster. 1V  Concentrator  control  should  provide  optimum 
energy  transfer  to  the  solar  thruster  while  protecting  the 
concentrator  and  spacecraft  from  the  dynamic  behavior  of 
the  whole  system. 

Concentrator  control  should  provide  optimum  energy  transfer 
to  the  solar  thruster  while  protecting  the  concentrator  and 
spacecraft  from  the  dynamic  behavior  of  the  whole  system. 

Focal  Spot: 

Proper  tracking  of  the  Sun  or  proper  positioning  of  the 
focal  spot  ensures  maximum  energy  transfer  to  the  fluid 
being  heated.  To  properly  track  the  Sun  for  maximum  heat 
transfer,  the  controller  requires  a  sensor  with  enough 
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resolution  to  ensure  proper  control  or  better  optimal 
control  for  maximum  heat  transfer. 

The  resulting  sensor  is  shaped  like  a  cone  and  serves 
not  only  as  the  sensor  but  also  as  the  solar  heat  absorber. 
The  method  of  determining  where  the  focal  spot  is  located 
using  the  new  sensor  system  is  based  on  the  method  used  in 
a  Shack-Hartmann  optical  wave  front  sensor.  The  solar 
sensor  developed  in  this  dissertation  could  be  applied  to 
either  ground  or  space-based  systems 


Chapter  2  Solar  Sensor  Options'7  Vl 

Most  terrestrial  Sun  trackers  were  developed  to  track 
the  Sun  from  the  Earth  and  rely  on  open  loop  mechanical 
tracking  algorithms.  That  is  tracking  then  was  based  on 
mechanical  position  and  vertical  position  error.  Relying  on 
open  loop  control,  these  systems  were  prone  to  error, 
because  they  lacked  Sun  position  error  and  inertial 
reference  feedback.  The  open  loop  control  systems  are  also 
unable  to  handle  the  strut  length  thermal  dependence  that 
feedback  systems  are  able  to  overcome. 

2.1  Thermal  based  sensors,  such  as  the  pyranometer  and 
pyrheliometer  rely  on  thermopiles  or  thermocouples  to 
measure  the  temperature  rise  or  difference  between  the 
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sensor  and  its  environment.  The  pyranometer  uses  an 
alternating  black  and  white  plate  connected  to  the 
thermopile.  The  incident  solar  radiation  is  reflected  and 
absorbed  differently  resulting  in  different  temperatures  on 
the  black  and  white  areas.  The  thermopile  then  can  measure 
the  temperature  differences  between  the  black  and  white 
areas  as  the  measurement  of  the  solar  radiation.  Another 
version  of  the  pyranometer  uses  a  Silicon  diode  as  the 
radiation  detector.  The  diode  produces  a  current  directly 
related  to  the  solar  radiation.  With  the  photodiode,  the 
pyranometer  is  inexpensive,  insensitive  to  temperature  or 
tilt  effects  and  does  not  reguire  an  external  power  supply. 

2.2  Normal  Incident  Pyrheliometer  (NIP):  These  devices 
are  similar  to  the  pyranometer  except  that  they  require  a 
Sun  electrically  driven  equatorial  mount  to  follow  the  Sun 
for  measuring  direct  sunlight.  The  unit  consists  of  a 
thermopile  and  a  collimator  tube  that  limits  incident 
radiation  to  a  5  degree  circular  cone.  Incident  radiation 
heats  the  thermopile  which  converts  the  temperature  rise 
into  a  voltage  representing  the  incident  radiation. 

2.3  Black  body  cavities  have  also  been  used  as 
absorbers,  storage  devices  and  heat  exchangers.  Black 
bodies  are  defined  in  most  physics  books  as  a  body  that 
absorbs  all  radiation  that  strikes  it.  Although,  no  true 
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black  body  exists  a  graphite  body  is  capable  of  absorbing 
all  but  3  percent  of  incident  radiation  striking  it.  Good 
approximations  of  black  bodies  are  cavities  having  a  small 
aperture  or  opening  through  which  the  radiation  enters  the 
cavity.  The  radiation  cannot  escape  out  of  the  box  and 
increases  the  box  temperature.  These  devices  measure 
radiation  with  either  a  photodiode  or  thermocouple  mounted 
in  the  cavity.  Once  the  temperature  of  the  black  body  is 
known,  the  radiation  can  be  determined  from  Stephan- 
Boltzmann' s  law. 

2.4  Photo  Diodes  and  Photo  Cells  are  electronic 
sensors  that  can  be  used  as  replacements  for  thermocouples 
or  thermopiles  in  other  devices.  Photo  diodes  consist  of  a 
Silicon  diode  in  a  container  that  exposes  the  PN  junction 
to  incident  light.  These  devices  operate  on  the  internal 
photoelectric  effect.  Photo  cells  are  also  Silicon-based 
and  are  usually  operated  at  short-circuit  current.  Both 
photo  cells  and  photo  diodes  are  used  to  detect  visible 
light . 

2.5  Net  Radiometers  are  devices  that  measure  all¬ 
wavelength  radiation  by  registering  the  temperature 
determined  by  the  balance  between  absorbed  and  radiated 
flux  on  two  horizontal  surfaces  thermally  separated  from 
each  other.  The  top  surface  of  the  radiometer  has  a 
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coating,  usually  white,  that  a  thermocouple  or  thermopile 
is  connected  along  with  a  temperature  measuring  device 
connected  to  the  bottom  surface  to  measure  the  two 
temperatures  to  determine  the  difference  between  reflected 
radiation  and  absorbed  radiation. 

2 . 6  Calorimeters  are  devices  that  measure  a 
temperature  increase  or  difference  to  indicate  radiant 
energy  from  the  Sun.  The  simplest  type  of  calorimeter  is 
the  "coffee  cup"  style  which  is  simply  two  Styrofoam  cups 
inside  each  other  with  water  and  a  thermometer  in  the  water 
to  measure  the  temperature  increase.  The  outside  cup 
serves  as  an  insulator  for  the  inner  cup  so  that  the 
temperature  increase  in  the  water  can  be  attributed  only  to 
the  energy  transferred  to  the  water.  The  medium  that  is 
between  the  two  cups  is  either  air  or  vacuum,  with  vacuum 
providing  a  better  insulator  similar  to  a  Dewar.  Other 
types  of  calorimeters  exist  and  use  alternatives  to  the 
water  for  the  mass  that  increases  in  temperature.  Some 
calorimeters  consist  of  a  metal  chamber  with  an  inner  metal 
chamber  for  the  water  or  thermal  mass.  A  temperature 
sensing  device  such  as  a  thermometer  passes  through  the 
inner  and  outer  chamber  to  register  the  temperature  change. 

Chapter  3  Solar  Sensor  Problem  Statement 


li 


To  be  useful  as  a  solar  sensor  both  on  the  Earth  and 


in  space  the  sensor  cannot  rely  on  Earth  coordinates  for 
tracking  information  and  needs  to  be  robust  enough  to 
provide  solar  tracking  while  the  thruster  is  maneuvering. 
The  sensor  also  must  be  able  to  handle  the  high 
temperatures  involved  with  a  solar  concentrator  that 
concentrates  the  sunlight  to  radiation  limited  temperatures 
approaching  4000  degrees  Kelvin  on  the  thruster. 

The  research  problem  is  to:  (1)  develop  a  system  that 
provides  information  about  the  location  of  the  focus  beam 
from  the  concentrator  that  also  can  tolerate  the  extreme 
temperatures  involved  in  providing  solar  heating,  and  (2) 
develop  the  mathematics  of  a  control  system  that  will  focus 
the  concentrator  in  an  optimal  fashion  based  on  the  sensor 
measurements.  The  system  presented  in  this  dissertation 
develops  a  new  sensor  that  accomplishes  both  of  the 
requirements  for  a  ground  based  and  a  space  based  solar 
control  system. 

NIPs  cannot  be  used  for  a  thruster  sensor  as  they 
track  the  Sun  mechanically  and  are  tied  to  Earth 
coordinates  for  locating  the  Sun.  NIPs  as  well  as 
pyranometers  cannot  function  at  the  temperatures  developed 
by  the  concentrator.  Also,  photodiodes  and  photo  cells  do 
not  have  the  temperature  capabilities  necessary  to  directly 
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be  used  as  a  focal  spot  sensor.  Photo  devices  are  limited 
to  a  temperature  of  about  85  degrees  C  while  the 
temperature  capability  of  the  pyranometer  is  about  40 
degrees  C. 

A  black  body  or  cavity  radiation  sensor  would  not  be 
useful  to  sense  the  focal  spot  as  the  thruster  does  not 
have  a  cavity  on  board.  Also,  the  black  body  or  cavity 
suffers  from  a  lack  of  response  time,  being  too  slow  for 
direct  thrust.  The  black  body  or  cavity  sensor  is  better 
suited  to  a  thruster  that  utilizes  stored  energy  for 
propulsion  and  does  not  have  a  need  for  continuous  thrust. 

Calorimeters  suffer  many  of  the  same  problems  that  the 
other  sensors  suffer.  Calorimeters  need  to  be  in  the 
direct  solar  radiation  path  to  function  so  would  in  fact 
impair  or  block  the  sunlight  from  reaching  the  thruster. 
Also,  the  materials  used  in  constructing  the  calorimeter 
would  limit  the  temperatures  that  it  could  withstand  to 
several  hundred  degrees  C  instead  of  the  required  3100 
degrees  K  required  in  a  concentrated  solar  light  situation. 
The  need  for  a  thermal  mass  represented  by  the  water  in  the 
coffee  cup  calorimeter  limits  the  feasibility  of  using  such 
a  device  on  a  spacecraft.  Calorimeters  can  not  be  used  in 
a  continuous  temperature  update  mode  that  limits  them  from 
use  as  a  focal  spot  sensor.  The  resolution  and 
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directionality  required  from  the  sensor  is  another 


limitation  of  the  calorimeter  as  they  cannot  provide 
directional  information  needed  for  focal  spot  tracking. 

Radiometers  are  another  type  of  device  that  have 
temperature  limitations  that  would  disqualify  them  from 
being  used  as  focal  spot  sensors.  The  maximum  temperature 
that  can  be  withstood  by  radiometers  is  40-85  degrees  C. 

The  sensor  system  to  be  developed  solves  the 
temperature  issue  by  using  the  absorber  as  a  wave  front 
sensor.  Since  the  absorber  has  to  take  the  high 
temperatures  of  the  concentrated  sunlight,  it  is  already 
capable  of  handling  the  temperature  requirements.  The  CCD 
images  the  absorber  to  determine  the  location  of  the  focal 
beam  and  may  use  thrust  as  an  additional  parameter  to 
improve  the  focal  beam  tracking  performance. 


Chapter  4  Experimental  Research 

4.1  Fourier  Analysis  and  Focal  Spot  Simulation 

Experiment  1 

The  appearance  of  the  distribution  of  intensity  in  the 
focal  spot  is  modified  by  the  surface  it  is  scattered  from, 
i.e.,  the  absorber.  The  detected  intensity  is  complicated 
by  shadows  and  the  specular  component  of  reflection.  It  is 
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easy  to  recognize  the  absorber  but  difficult  to  determine 
where  on  the  absorber  most  of  the  energy  is  going  using 
just  human  perception.  Furthermore,  at  different  points  in 
the  focal  spot,  the  direction  of  incident  light  varies. 
Therefore,  as  the  focal  spot  drifts  across  the  absorber  the 
distribution  of  scattered  light  changes  in  a  complicated 
way  and  it  is  difficult  to  determine  how  to  correct 
pointing.  This  difficulty  has  been  demonstrated  in  the 
laboratory  in  concentrator  testing. 

Digital  signal  processing  can  be  used  to  "de-convolve" 
the  misaligned  focal  spot  image  from  the  scattering 
distribution  of  the  thruster  absorber.  In  the  case  of  the 
current  section,  the  sensors  used  for  focal  spot  placement 
consist  of  a  coarse  alignment  sensor  for  gross  movement  of 
the  collectors,  followed  by  fine  alignment  sensors  for 
precise  alignment  of  the  focal  spot.  Proper  alignment  is 
determined  by  examining  the  shape,  size  and  intensity  of 
the  light  of  the  focal  spot  on  the  absorber. 

Sensors  that  allow  this  type  of  examination  and  that 
would  handle  the  high  temperatures  presented  at  the 
absorber  are  not  readily  available.  Two  possibilities  are 
Charge  Coupled  Device  (CCD)  arrays  and  photo  diode  detector 
arrays.  The  CCD  array  would  be  positioned  out  of  direct 
sunlight,  using  a  lens  system  for  imaging,  and  would 
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therefore  be  somewhat  protected  from  the  intense  heat.  The 
photo  diodes  as  discussed  in  section  3  would  be  in  direct 
sunlight  and  would  not  be  able  to  handle  the  high 
temperatures . 

To  begin  the  work,  a  computer  simulation  of  the  focal 
beam  of  an  off-axis  concentrator  was  used  to  explore  image 
methods  to  use  on  the  images  for  determining  the  focal  spot 
location . 

Computer  Simulation: 

A  C  program  called  Off-axis  generates  the  focal  spot 
information,  based  upon  having  a  Lambertian  target  at  the 
focal  plane.  This  program  was  written  by  Dr.  Michael  R. 
Holmes  of  the  Air  Force  Research  Laboratory,  specifically 
for  inflatable  concentrators,  and  provides  an  "image"  of 
the  focal  spot  target  area  being  illuminated  by  an  off-axis 
concentrator.  Output  data  generated  by  offaxis.c 
represented  the  focal  spot  image  intensity. 

Although,  the  simulation  produces  ideal  focal  beams, 
using  ray-tracing  with  Monte  Carlo  slope  error  simulation, 
it  is  an  important  starting  place  to  prove  that  the 
concepts  in  this  paper  could  work,  and  supplied  important 
data  for  analysis  to  begin. 

Since  the  simulation  produced  a  nearly  ideal 
situation,  the  focal  spot  generated  resulted  in  a  Gaussian- 
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type  focal  spot  intensity  distribution.  Contour  and  surface 
plots  of  the  simulation  results  are  shown  in  Figures  4-7 
indicating  the  Gaussian  3-dimensional  nature  of  the  focal 
beam. 

3-D  Plot  Offset  Angle  Jhat  =  0.5 


Figure  4  3D  PLOT  OF  FOCAL  SPOT  JHAT  -0.5 


Figure  5  CONTOUR  PLOT  OF  FIGURE  2 
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3-D  Plot  Offset  Angle  ihat  =  0.5 


Figure  6  3D  PLOT  OF  FOCAL  SPOT  IHAT  =0.5 


Contour  Plot  Offset  Angle  ihat  =0.5 
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Figure  7  CONTOUR  PLOT  FOR  FIGURE  4 

The  contour  plots  show  that  the  focus  spot  is  actually 
elliptical  instead  of  perfectly  circular.  The  elliptical 
nature  of  the  focus  is  due  to  the  off  axis  nature  of  the 
concentrator  and  will  be  utilized  in  the  determination  of 
the  location  of  the  focus.  The  parameters  for  the  Sun  and 
concentrator  positions  are  passed  to  the  program  through 
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of f axis . h . 


By  varying  the  parameters  in  this  file,  the 


various  misaligned  positions  of  the  Sun  and  concentrator 
can  be  simulated.  Of  primary  concern  to  this  section  are 
the  Sun's  offset  angle  parameters.  These  two  parameters 
change  the  angular  position  of  the  solar  disc  in  the  x, 
(IHAT)  and  y,  (JHAT)  directions.  See  Figure  8  for  the 
geometry . 

Solar  Thruster  Concentrator; 


Figure  8  OFF-AXIS  PARABOLOID  GEOMETRY 


The  simulation  runs  for  this  portion  of  the 
dissertation  varied  these  two  parameters  only.  Note  that 
once  the  data  was  generated,  the  image  matrix  was  treated 
in  an  X-Y  image  data  plane  instead  of  what  is  represented 
in  Figure  8  that  is  the  focal  spot  data  image  is  presented 
in  its  own  X-Y  plane  separate  from  the  coordinate  system 
depicted  in  Figure  8  and  depicted  in  Figure  9.  The 
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specific  concentrator  modeled  in  the  runs  was  a  4  X  6  meter 
Flight  Scale  Concentrator  manufactured  by  SRS  Technologies, 
Huntsville,  AL . 


Vector  From  CCD  or 
T  a  r get  to  Pixel  on 
Concentrator  under 
Consideration 


Vector  from  origin 
to  CCD  location 


Image  Coordinates 
u  =  AuTe*, 
v  =  Av*j’(T2 


Figure  9  CCD  or  Imager  Coordinates 

Analysis  Methods: 

Once  the  image  data  were  generated,  the  data  were 
analyzed  using  Matlab.  Initial  analysis  used  the  2- 
dimensional  Fourier  Transform  (FT)  or  Fast  Fourier 
Transform  (FFT)  as  Matlab  implements  this  transform. 
Equation  1  shows  the  definition  of  the  2-dimensional 
continuous  Fourier  Transform,  with  the  limits  of 
integration  being  ±o° .  This  form  of  analysis  was  used  as  a 
starting  point  as  it  was  based  upon  previous  work  in  one 
dimensional  signal  processing  and  no  a  priori  expectations 
about  what  would  work.  The  one  dimensional  signal 
processing  used  Fourier  techniques  that  are  expanded  to  two 
dimensional  or  images  quite  readily. 
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f  Jf(x,yyj{aK+^dxdy 


The  FFT  is  a  fast  mathematical  algorithm  for  computing 
the  discrete  form  of  Equation  1;  see  Equation  2  for 
discrete  FT,  utilized  by  computers  to  speed  up  the 
calculation  of  FT's  in  one  dimension.  In  general  cases, 
the  2-dimensional  FFT  may  be  calculated  by  using  the  1- 
dimensional  case  and  applying  it  to  the  data  twice. 

X{k)  =  Yn--JWe~imNnk  (2) 

with  coN  =  2jt/N,  and  summed  over  n. 

However,  this  effort  turned  out  to  be  futile  as  the 
Fourier  Transform  of  a  Gaussian  shaped  signal  is  also 
Gaussian.  See  Figure  10. 

Plot  of  Focal  Beam  offset  angles  =  0.0 
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Figure  10  3D  PLOT  OF  FFT  FOCAL  SPOT  BEAM 
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Although,  in  the  ideal  modeling  world  represented  by 
the  focal  beam  data  from  the  simulation  a  maximum  location 
could  be  found,  i.e.,  the  beam  maximum  found  using  a 
gradient  operation  on  the  original  data.  A  method  that 
could  handle  a  visually  more  complex  image  with  specular 
reflections  co-mingled  with  a  non-ideal,  Gaussian  shaped 
focal  spot  would  need  to  be  developed  to  be  useful  in  a 
spacecraft  or  real  case. 

The  next  method  of  analysis  performed  was  a  variant  of 
the  Short  Time  Fourier  Transform  (STFT)  for  time  signals. 

The  image  data  was  treated  as  a  matrix  in  Matlab  and  a 
vector  Fourier  Transform  (FT)  was  performed  on  the  matrix. 
Matlab  performs  the  FFT  on  a  matrix  by  performing  the 
transform  on  each  column.  What  this  does  to  the  matrix,  in 
effect,  is  perform  a  windowed  FFT  or  STFT  on  the  image. 

See  Equation  3,  w(t)  is  the  window  function  in  the  1 
dimensional  continuous  transform. 

oo 

F((d)  =  Jw(t  -r)f(t)e~jco,dt  (3) 

—  00 

In  fact,  the  window  happens  to  be  rectangular  shaped  and  is 
exactly  one  column  size  in  length.  The  shift  value, x  for 
each  FFT  also  happened  to  be  one  column  length,  thus  making 
the  calculation  fairly  straightforward  to  accomplish  and 
actually  located  at  each  X  value.  The  result  of  using  the 
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STFT  on  the  image  was  a  2-  dimensional  graph  depicting  the 
FT  of  the  image  in  lines  along  the  Y  axis  on  grid  spacing 
identical  to  the  X  axis  grid.  Thus,  the  STFT  provided  a 
freguency  distribution  along  the  Y  axis,  but  local  to  each 
X  value. 

The  position  of  the  focal  beam  was  determined  by 
finding  the  maximum  value  of  the  transforms  in  each 
direction.  These  maximums  were  then  used  to  calculate  the 
amount  of  correction  in  each  direction  needed  to  re-center 
the  beam.  The  FFT  and  STFT  approach  was  chosen  to  see  what 
could  done  with  the  images  with  no  a  priori  expectations. 

Results  of  Experiment  1: 

The  results  of  performing  the  STFT  on  various  offset 
angles  are  shown  in  Figures  11-16  .  Comparing  the  data 
plots  from  the  STFT's  to  the  contour  plots.  Figures  17-19 
of  the  same  offset  angles,  gave  the  sense  that  STFT'' s  would 
provide  the  needed  information  about  the  location  of  the 
focal  spot  in  the  X  and  Y  directions.  Table  1  describes 
each  of  the  graphs . 


Figure  11 

STFT  X  Dir 

I HAT  =  -0 . 5deg. 

Figure  12 

STFT  Y  Dir 

I HAT  =  -0 . 5deg. 

Figure  13 

STFT  X  Dir 

JHAT  =0.5  deg. 

Figure  14 

STFT  Y  Dir 

JHAT  =0.5  deg. 

Figure  15 

STFT  X  Dir 

IHAT=JHAT=0 . 1 

deg. 
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Figure  16 

STFT  Y  Dir 

IHAT=JHAT=0 . 1 

deg. 

Figure  17 

Contour 

IHAT  =  0.5  deg. 

plot 

Figure  18 

Contour 

JHAT  =0.5  deg. 

plot 

Figure  19 

Contour 

IHAT=JHAT=0 . 1 

plot 

deg. 

Table  1 


3-D  Plot  STFT  Offset  Angle  ihat  =  0.5 


Figure  11  STFT  X  DIRECTION,  IHAT  -0.5 
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Figure  12  STFT  Y  DIRECTION,  I HAT  =  -0.5 


3-D  Plot  STFT  X  Dir.  jhat  =  0.5 
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Figure  13  STFT  X  DIRECTION,  JHAT  =0.5 
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3-D  Plot  ST  FT  Y  Dir.  jhat  =  0.5 


Figure  14  STFT  Y  DIRECTION,  JHAT  =0.5 

STFT  X  Direction  Offset  Angle  jhat  =  0.1 


Figure  15 


STFT  X  DIRECTION, 


JHAT 
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STFT  X  Direction  Offset  Angle  jhat  =  0.1 


Figure  16  STFT  Y  DIRECTION,  JHAT  =  I HAT  =0.1 


Contour  Plot  Offset  Angle  ihat  =0.5 


X  Direction  Pixel  Number 


Figure  17  CONTOUR  PLOT  IHAT  =0.5 
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Contour  Plot  Offset  Angle  jhat  =  0.5 


X  Direction  Pixel  Number 


Figure  18  CONTOUR  PLOT  JHAT  =0.5 


Contour  Plot  Offset  Angles  =  0.0 


Figure  19  CONTOUR  PLOT  JHAT  =  I HAT  =0.1 

Computing  the  STFT  first  along  the  Y  direction 

(columns  of  the  matrix)  determined  position  of  the  beam  in 
the  X  direction.  Computing  the  STFT  along  the  X  direction 
(using  the  transposed  matrix)  determined  the  location  of 
the  focal  beam  in  the  Y  direction.  Finding  the  maximum  of 
the  spectrum  and  correlating  that  to  the  pixel  number  of 
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that  maximum  gave  the  X,  Y  location  of  the  focal  spot  of 
the  focal  beam.  Table  two  gives  the  results  of  finding 
maximum  values  for  each  of  the  cases  in  Figures  11-16  . 


Figures  11- 

12 

X  pixel 

Number  =  76 

Y  pixel 

Number  =  90 

Figures  13- 

X  pixel 

Y  pixel 

14 

Number  =  56 

Number  =  76 

Figures  15- 

X  pixel 

Y  pixel 

16 

Number  =  72 

Number  =  76 

Table  2 


Each  of  the  pixel  numbers  in  table  2  would  be 
converted  to  physical  measurements  using  the  scaling 
factors  in  of f axis. c.  The  factor  in  the  computer  runs 
presented  in  this  dissertation  is  0.226  units  of 
measurement  per  pixel.  The  units  in  most  cases  are  either 
centimeters  or  millimeters. 

A  problem  with  the  use  of  the  STFT  in  the  situations 
presented,  is  that  some  of  the  offset  angles  produced  only 
a  change  in  magnitude  and  the  width  of  the  focal  beam  and 
not  necessarily  a  change  in  its  X,  Y  location.  The  beam 
was  not  expected  to  change  peak  frequency  value,  but  it  was 
expected  to  change  position  in  either  X  or  Y.  When  the 
beam  expanded  or  narrowed  due  to  the  misalignment  but  did 
not  move  the  beam  in  X  or  Y,  the  STFT  did  not  provide  the 
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necessary  position  information.  The  STFT,  however,  was 


able  to  show  that  the  focal  beam  width  and  the  maximum 
value  of  the  focal  beam  had  changed,  but  only  when  the 
offset  plot  was  compared  to  the  zero  offset  plots, 
visually . 


Summary  of  Experiment  1 : 

Data  from  a  computer  program  was  used  to  generate 
nearly  ideal  concentrator  performance  that  was  analyzed 
using  the  STFT.  The  results  of  using  the  1  dimensional 
STFT  on  images  showed  promise,  thus  verifying  the  concepts, 
in  using  the  STFT  in  determining  the  placement  of  the  focal 
beam.  The  STFT's  were  performed  first  along  the  Y  axis  to 
determine  the  pixel  number  of  the  maximum  beam  value  in  the 
X  direction,  then  the  STFT' s  were  performed  along  the  X 
axis  to  determine  the  pixel  number  of  the  maximum  beam 
value  in  the  Y  direction.  These  results  provided  incentive 
to  continue  analyzing  the  data  using  the  STFT. 

One  problem  that  occurred  when  using  the  STFT  was  that 
in  some  instances,  the  beam  only  changed  in  magnitude  and 
not  direction.  In  this  case,  the  STFT  does  not  provide 
sufficient  information  to  get  control  information. 

Using  the  STFT  on  simulated  data  representing  the 
intensity  distribution  of  the  focal  point  was  instrumental 
in  showing  that  the  concept  of  imaging  the  absorber  with 
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light  on  it  could  be  used  to  determine  the  position  of  the 


spot . 

Once  the  location  of  the  focal  spot  is  determined, 
control  inputs  should  be  calculated  for  output  to  the 
hexapod  to  re-position  the  spot  for  better  performance. 

The  next  section  builds  on  the  present  section'’  s  results  by 
imaging  an  actual  thruster  and  applying  the  STFT  to  the 
images  along  with  some  research  into  image  enhancements 
that  may  apply  to  the  images  used  for  control. 


4.2  Image  Processing  of  Thrusters: 

Since  control  over  the  focal  spot  location  on  the 
thruster  absorber  is  paramount  in  delivering  maximum  heat 
to  the  propellant,  it  drives  the  search  for  a  method  of 
locating  the  focal  spot  on  the  thruster  and  deriving 
concentrator  control  inputs  to  reposition  the  focal  spot  to 
optimize  power  or  heat  transfer  to  the  propellant.  The 
sensor  and  data  analysis  algorithms  used  to  determine  focal 
spot  location  are  the  most  important  parts  in  the  overall 
control  mechanism.  The  sensor  used  in  the  experiment 
performed  in  this  section  of  this  chapter  was  a  CCD  camera 
imaging  a  thruster. 
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Previous  work  on  this  topicvi;L  and  presented  in  the 
previous  section,  indicated  that  a  short  time  Fourier 
transform  (STFT)  on  the  image  data  would  lead  to  the 
required  focal  spot  location  information  for  data  from  a 
computer  simulated  off-axis  concentrator.  Also, 
information  in  that  section  indicated  that  the  STFT  was  not 
sufficient  in  cases  where  the  movement  of  the  concentrator 
did  not  move  the  focal  spot  in  X  and  Y,  but  only  changed 
the  maximum  value.  The  purpose  of  this  effort  is  to  move 
from  the  computer  simulation  to  actual  thruster  datavli:L. 

Experiment  2: 

The  setup  for  the  concentrator  and  data  collection 
system  is  shown  in  Figure  20.  A  1  m  X  2  m  elliptical 
concentrator  was  utilized  for  the  concentrator  as  that  was 
the  style  of  concentrator  available  at  SRS  Technologies  in 
Huntsville,  AL .  Figures  85  and  86  in  appendix  A  show 
photos  of  the  actual  concentrator  and  light  source.  This 
type  of  concentrator  was  acceptable  as  the  data  that  was  to 
be  collected  was  only  CCD  camera  shots  of  the  thruster  on 
"Sun"  and  only  a  focal  spot  on  the  thruster  was  needed. 

The  experimental  setup  in  Figure  20  and  discussed  in 
this  section  was  optically  equivalent  to  the  situation 
encountered  in  space  utilizing  larger  concentrators.  A 
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divergent  point  source  located  at  the  other  focus  of  the 
ellipse  was  used  to  simulate  the  light  from  the  Sun. 

Because  the  concentrator  was  elliptic,  the  divergent  source 
could  be  located  at  the  far  focus  (because  light  rays 
emitted  from  one  focus  of  an  ellipse  get  reflected  through 
the  other  focus)  while  the  thruster  was  located  at  the  near 
focus  as  shown  in  the  Figure. 

The  setup  in  Figure  20  allowed  experiments  to  be 
performed  in  a  controlled  environment  that  allowed  various 
techniques  to  be  tried  before  used  on  a  spacecraft.  Figure 
21  shows  a  schematic  (2  Dimensional)  of  the  setup  in  Figure 
20.  Appendix  A  presents  photographs  of  all  of  the 
equipment  and  setup.  A  scissor  jack  platform  was  used  to 
vary  the  vertical  position  of  the  point  source  at  the  far 
focus.  Figure  21  also  shows  the  axes  positioning  for  the 
source,  which  is  then  used  in  the  production  of  the  images 
presented  later.  The  zero  height  position  of  the  source  is 
at  the  focus  location  of  16.25  inches  above  the  floor  of 
the  setup  building. 

The  data  collection  system  consisted  of  an  ST-6  CCD 
Camera  manufactured  by  SBIG  in  Santa  Barbara,  CA  and  a  Sony 
Vaio  computer  running  SBIG' s  proprietary  image  collection 
software.  This  software  allowed  the  user  to  control  and 
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set  various  parameters  of  the  camera.  One  parameter  under 
user  control  was  the  exposure  time  for  the  camera  to  view 
the  scene.  Exposure  times  were  important  as  they 
determined  whether  the  signal-noise  ratio  of  the  image  was 
sufficient  to  provide  a  good  image  and  not  saturate  the 
image.  The  images  presented  later  in  this  paper  indicate 
the  exposure  time. 

The  camera  used  a  40  mm  f.l  lens  with  a  small  aperture 
similar  to  a  pinhole  to  provide  good  depth  of  field.  The 
CCD  was  spaced  from  the  lens  to  give  the  best  full  view 
image  of  the  thruster  in  focus  at  1  meter  from  the 
thruster.  The  CCD  in  the  camera  measured  23  X  27  micro¬ 
meters  and  had  375  pixels  by  241  pixels.  As  Figure  20 
presents,  the  camera  was  mounted  on  the  concentrator 
supporting  ring  approximately  64  inches  away  from  the 
thruster  surface. 
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Test  Apparatus 


Figure  20  TEST  SETUP 


Schenatic  (ncrt  to  scole) 


Figure  21  2  D  SCHEMATIC  OF  TEST  SETUP 

Method  of  Data  Collection: 

Since  the  algorithms  must  be  able  to  resolve  changes 
in  translation  of  0.1  inch  and  0.1  degrees,  the  calculation 
for  an  angle  of  0.1  degree  was  performed  for  the  setup 
shown  in  Figure  20,  and  was  found  to  be  about  1  inch. 
Therefore,  the  source  light  was  moved  in  1  inch  increments 
in  both  the  vertical  and  horizontal  directions.  A  standard 
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measuring  tape  was  used  to  measure  the  placement  of  the 
source  light;  measurements  were  measured  to  1  inch  lines 
and  could  be  +  or  -  an  eighth  of  an  inch.  Once  the 
position  of  the  source  was  set,  the  camera  was  commanded  to 
take  an  image  of  the  thruster.  Each  image  included  taking 
a  "dark"  image  with  the  camera  element  blocked  to  set  a 
zero  image  that  is  subtracted  from  each  light  image.  Each 
image  was  taken  for  0.1  seconds.  At  the  time  that  the 
experiment  was  performed,  it  was  thought  that  the  analysis 
to  be  performed  would  provide  enough  information  regarding 
the  specular  reflections,  that  no  other  exposure  times  were 
utilized.  Table  3  shows  each  of  the  data  runs.  The  maximum 
vertical  limit  of  the  scissors  jack  was  4  5/8  inches  while 
the  maximum  horizontal  position  was  3  inches  off  of  focus. 
The  horizontal  limit  was  a  result  of  how  much  the  scissors 
jack  and  point  source  could  be  moved  on  the  block  before 
the  unit  would  fall  off  of  the  block. 


Run  Number 

Offset  Vertical 

(inches) 

Offset  Horizontal 

(inches) 

Exposure  Time 

Comments 

1 

0 

0 

0 . 1  s 

Start 

2 

1 

0 

0.1s 

3 

2 

0 

0.1s 
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Run  Number 

Offset  Vertical 

Offset  Horizontal 

Exposure  Time 

Comments 

(inches) 

(inches) 

23 

4 

3 

0 . 1  s 

Table  3:  Test  information  for  experiment. 


Image  Processing: 

Mathworks  MATLAB  software  with  the  Image  Processing 
Toolbox  was  used  to  perform  the  image  analysis.  Each  image 
was  loaded  into  the  MATLAB  program  along  with  the  one 
source-off  image  taken.  The  source-off  image  was  taken  to 
capture  the  ambient  lighting  effects  in  the  SRS  high  bay. 

The  methods  used  to  process  the  images  consisted  of 
several  steps.  The  methods  used  and  described  next  were 
taken  from  those  methods  discussed  in  computer  vision 
texts.  Many  computer  vision  texts  develop  methods  that  try 
to  mimic  the  human  eye  and  mimic  human  interpretation  of 
those  images.  Thus  histogram  equalization  was  thought  to 
be  a  necessary  enhancement  to  perform  on  the  images  as  that 
would  improve  a  visual  interpretation  of  those  images. 

Also,  it  was  thought  that  it  would  be  necessary  to 
determine  or  locate  the  edges  of  each  of  the  tubes  to 
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facilitate  locating  the  images  of  the  concentrator  on  the 


sensor . 

Each  image  was  first  enhanced  to  provide  contrast 
adjustment  so  that  a  human  could  interpret  the  image  and 
locate  various  portions  of  the  sensor.  Contrast 
enhancement  was  thought  to  be  necessary  as  each  image  had  a 
very  limited  range  of  intensities,  as  evidenced  in  a 
representative  histogram  shown  in  Figure  23,  so  the  images' 
ranges  of  intensities  were  stretched  to  fill  the  dynamic 
range.  Then  each  image  was  filtered  using  a  Gaussian  low- 
pass  filter.  After  low-pass  filtering,  each  image  was  then 
filtered  using  a  LaPlacian  high-pass  filter  (LOG— LaPlacian 
of  Gaussian)  to  provide  image  edge  enhancement.  The  LoG 
filtering  was  thought  to  help  eliminate  the  histogram 
problems  encountered  in  the  images  taken  for  the 
experiment . 

Although,  the  LoG  filter  removes  intensity  bias, 
caused  by  changes  in  illuminations,  it  did  not  help  improve 
the  histogram  of  the  imageslx.  This  last  image,  in  effect  a 
band-passed  image,  with  the  edges  enhanced  and  the  noise 
from  the  LaPlacian  filter  reduced  by  the  low-pass  filtering 
was  used  for  further  analysis. x  Finally,  the  pre-processed 
and  filtered  source-off  image  was  subtracted  from  each 
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image  to  remove  the  effects  of  the  ambient  lighting  which 
could  add  to  any  bias  noted  in  the  image. 

The  second  method  used  the  same  pre-processing  steps 
as  the  first  one  but  substituted  an  average  filter  instead 
of  the  LOG  filter.  Each  method  will  be  discussed  as  to  its 
effectiveness  in  the  next  section. 

Image  Analysis: 

Figure  22  shows  an  image  with  the  source  located  at  a 
4-inch  vertical  offset  and  before  histogram  equalization. 

A  histogram  of  the  image  with  a  4-inch  vertical  offset  is 
shown  in  Figure  23  and  possibly  indicates  the  need  to 
perform  contrast  enhancement,  as  most  of  the  pixels  lie  in 
the  lower  portion  of  the  histogram.  The  need  for  contrast 
enhancement  was  in  line  with  the  thought  of  using  elements 
and  theories  from  computer  vision  to  locate  and  track  the 
focal  spot.  As  will  be  shown  later  in  this  dissertation, 
the  most  important  pixel  brightness  levels  are  located  at 
the  far  right  of  the  histogram  images,  representing  the 
bright  white  specular  reflections  of  the  concentrator. 
Figures  24  and  25  show  the  image  and  histogram  after 
histogram  equalization.  Figure  26  shows  the  source-off 
image  after  enhancement.  Figure  27  shows  the  final  image 
after  LOG  filtering  and  subtraction.  It  is  this  final 
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image  that  will  be  used  in  the  image  analysis.  Figure  28 
shows  the  final  enhanced  image  when  the  average  filter  is 
used  for  the  enhancement  process.  Each  of  the  images  shown 
in  Figures  27  and  28  were  then  analyzed  with  the  STFT  to 
find  out  whether  the  location  of  the  focal  spot  could  be 
determined.  Figure  29  shows  the  result  from  using  the  STFT 
on  the  data  from  LoG  filtering  of  the  image. 

Results  of  Experiment  2: 

Since  LoG  filtering  is  basically  a  band  pass  operation 
on  the  images,  the  STFT  produces  the  band  pass  response  of 
the  image  and  does  not  provide  information  for  determining 
focal  spot  location  using  the  STFT.  The  STFT  results  from 
using  the  averaging  filter  were  more  interesting  as  far  as 
their  use  in  determining  the  location  of  the  focal  spot. 

In  both  the  X  direction  and  the  Y  direction,  pixel 
numbers  for  the  maximum  of  the  focal  spot  were  determined 
with  ease.  The  STFT  in  X  and  Y  directions  are  shown  in 
Figures  30  and  31  and  show  that  the  maximums  in  each 
direction  on  Figure  30  and  31  would,  respectively,  produce 
the  X  pixel  number  and  Y  pixel  number  of  the  location  of 
the  focal  spot.  The  information  provided  from  those  two 
Figures  calculations  could  then  be  processed  to  provide 
concentrator  control  commands . 
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A  conclusion  from  this  section  is  that  the  averaging 
filter  was  shown  to  be  better  suited  when  the  STFT  was  used 
to  determine  focal  spot  location.  Although  the  LoG 
filtering  was  not  as  useful,  it  could  still  be  useful  for 
procedures  that  do  not  perform  Fourier-type  operations, 
such  as  spatial-domain  operations  involving  gradients, 
feature  extraction,  etc. 

Histogram  processing  is  useful  when  each  image  is  to 
be  interpreted  by  a  human  to  make  sure  that  the  intensities 
fully  populate  the  histogram  thus  increasing  the  ability  of 
the  human  to  discern  where  the  tubes  are  in  the  picture  to 
help  move  the  focal  spot.  Instead  of  taking  images  at 
multiple  1-inch  increments,  fewer  increments  and  more 
exposures  should  be  taken  to  be  sure  that  the  histograms 
are  more  fully  populated.  A  study  of  the  specular 
reflection  seen  at  the  absorber  needs  to  be  accomplished, 
as  the  current  imaging  techniques  and  analysis  do  not 
properly  account  for  the  specular  reflections.  Plus 
filtering  techniques  beyond  LoG  filtering  would  need  to  be 
studied  to  determine  if  additional  filtering  techniques  are 
able  to  handle  the  specular  reflections  and  also  provide  a 
better  means  of  determining  the  location  of  the  focal  spot. 
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Image  4  inch  offset  (V)  before  equalization. 
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Figure  22  Image  Before  Equalization 

Figure  19A  Image  in  Figure  19  Before  Equalization 


Histogram  of  4  inch  offset  image 
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Figure  23  HISTOGRAM  4  INCH  VERTICAL  OFFSET 
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Image  4  inch  Offset  (V)  after  equalization 
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Figure  24  IMAGE  AFTER  HISTOGRAM  EQUALIZATION 


Histogram  4  inch  veritcal  offset  after  equalization 
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Figure  25  5X3  HISTOGRAM  AFTER  EQUALIZATION 


Figure  26  SOURCE-OFF  AFTER  EQUALIZATION 
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Final  image  after  filtering  and  subtraction 


Figure  27  FINAL  IMAGE  AFTER  LOG  FILTERING  AND  SUBTRACTION 


Final  Image  after  enhancement  and  source-off  image  is  subtracted. 
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Figure  28  FINAL  IMAGE  AFTER  ENHANCEMENT  AND  SUBTRACTION 


Y  STFT  mesh  plot 


Figure  29  PLOT  OF  STFT  Y  DIRECTION  LOG  FILTERED  IMAGE 
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X  STFT  mesh  plot 


Figure  30  X  DIRECTION  STFT  MAGNITUDE 

Y  STFT  mesh  plot 


Figure  31  Y  STFT  MAGNITUDE 

The  experiment  performed  and  reported  in  section  4.2, 
gave  promising  results  in  the  search  for  an  acceptable 
device  for  tracking  the  focal  spot  on  the  thruster. 

However,  another  approach  was  suggested  by  Dr.  Michael 
Holmes  that  still  used  CCD  camera  but  interpreted  the  image 
data  differently  than  before,  his  suggestion  also  did  not 
require  the  Fourier  transform  processing  of  the  images.  It 
also  did  not  necessarily  need  much  enhancement  of  the 
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resulting  images  as  discussed  in  section  4.2.  That  new 
approach  was  to  use  the  absorber  as  a  sensor  and  process 
the  images  similarly  to  a  Shack-Hartmann  wave  front  sensor. 
Section  4.3  introduces  the  new  sensor  and  algorithms. 

4.3  Wave  Front  Sensor Xl 

This  section  describes  how  the  absorber  is  to  be  used 
as  a  sensor  for  concentrator  control.  The  sensor  is  shaped 
like  a  cone  and  serves  not  only  as  the  sensor  but  also  as 
the  solar  heat  absorber.  The  method  of  determining  where 
the  focal  spot  is  located  on  the  sensor  uses  either  one  of 
two  methods  derived  from  a  Shack-Hartmann  optical  wave 
front  sensor. 

The  Hartmann  wave  front  sensor  was  developed  by 
Hartmann  in  1900  and  was  used  for  checking  optical 
telescopes  for  aberrations.  It  was  an  array  of  holes  in  a 
plate  placed  in  front  of  the  mirror  of  the  telescope  being 
checked  so  that  light  tubes  from  the  array  would  impinge  on 
the  telescope  mirror.  Two  photographic  plates  were  used  to 
collect  the  light  information.  One  plate  was  placed  just 
before  the  focal  spot  of  the  telescope  and  the  other  one 
was  placed  just  behind  the  focal  spot.  The  distance 
between  the  plates  was  chosen  so  that  the  light  rays  would 
be  separated  from  each  other.  The  path  of  the  light  rays 
was  traceable  by  measuring  the  centroids  of  the  two  images. 
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Thus,  Hartmann  was  able  to  determine  Figures  of  merit  for 
various  large  telescopes  . X11 

Dr.  Roland  Shack  developed  the  method  of  replacing  the 
holes  of  the  Hartmann  sensor  with  small  lenses  or  lenslets 
to  improve  the  capability  of  the  Hartmann  sensor  when 
taking  images  of  satellites  or  stars  at  the  same  time  as 
the  wave  front  phase  error  was  being  determined.  The  method 
was  developed  while  solving  a  problem  with  imaging 
satellites  and  stars  from  the  Earth.  Dr.  Aden  Mienel' s 
solution  was  to  determine  the  Optical  Transfer  Function  of 
the  atmosphere  at  the  same  time  as  the  image  being  taken. 
The  Hartmann  sensor  was  not  satisfactory  as  Mienel  could 
not  allow  the  array  to  cover  the  aperture  of  the  telescope 
because  it  would  eliminate  the  image  of  the  satellite  in 
favor  of  determining  wave  front  tilt  over  the  imaging  of 
the  satellites.  Mienel  tried  to  use  a  beam  splitter  to 
take  Hartmann  images  while  taking  images  of  the  satellites 
and  stars.  This  method  suffered  from  the  weak  intensity 
images  of  the  light  rays  in  the  Hartmann  sensor. 

Shack  suggested  that  replacing  the  apertures  in  the 
Hartmann  sensor  with  lenses  would  eliminate  the  problems 
that  Mienel  was  having  with  his  Hartmann  device.  The 
lenslets  focused  spots  onto  an  image  plane  such  that  the 
displacement  of  the  spots  from  the  ideal  position  of  the 
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Thus,  the  Shack- 


spots  indicated  wave  front  tilt.*111 
Hartmann  system  overcame  the  shortcomings  of  the  Hartmann 
sensor  while  increasing  the  sensitivity  of  the  Hartmann 
device  alone.  A  comparison  of  the  two  systems  is  depicted 
in  Figure  32.  The  Shack-Hartmann  sensor  is  not  suitable 
for  use  in  tracking  the  focal  spot  of  concentrated  solar 
rays  because  of  the  high  temperatures  encountered  in  a 
concentrated  sunlight  application. 


CCD 

Hartmann 

Sensor 


Front  of  plate 

Matrix  of  Apertures  or  Holes 


Sub-aperture  rocal  spots. 


Shack-Hartmann  ccd  Matrix  of  Lenslets 

Sensor 


Figure  32  SHACK-HARTMANN  VS.  HARTMANN 

Cylindrical  Mirrors 

Since  the  sensor  for  this  application  uses  coiled 
tubing  as  its  primary  configuration,  a  brief  review  of 
cylindrical  mirrors  is  necessary  to  understand  the 
reflections  generated  by  the  conical  tubing  absorber. 
Figure  33  shows  the  normal  situation  encountered  with  the 
use  of  a  cylindrical  mirror.  In  general  a  cylindrical 
mirror  will  tend  to  compress  an  object's  reflection  along 
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the  length  of  the  mirror. 


An  object's  reflection  also 


tends  to  get  larger  as  the  object  moves  towards  the  mirror. 


Concentrator  image 


. 


Cylindrical  mirror  (one  tube) 


Flat  mirror  in  vertical  dir. 
Image 


Object 

Angle  of 


1/Do  +1/Di  =  1/f 

mag  =  -DilDo 


Convex  mirror  in  horizontal  dir. 

fie  focal  pont 

f  is  0.5*  Rc 


Figure  33  CYLINDRICAL  MIRROR 
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Figure  34  Camera  View  of  Cylindrical  Mirror 
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random  un-focused  position  of 
concentrator. 


Figure  35  View  From  Sensor  Back  to  Concentrator 


Figure  36  Construction  for  Flux-tube  Calculation 

Figure  33  also  shows  the  reflections  from  the 
cylindrical  mirror  along  two  of  its  axes.  In  one  axis,  in 
the  straight  on  view,  the  cylindrical  mirror  acts  like  a 
flat  mirror  with  the  light  rays  adhering  to  the  rules  for  a 
flat  mirror.  That  is,  incident  light  is  reflected  such 
that  the  reflected  ray  returns  from  the  mirror  as  if  it 
originated  at  a  point  behind  the  mirror,  at  equal  distances 
to  the  object;  this  is  also  a  virtual  image.  Or  the 
reflection  occurs  such  that  the  reflected  ray' s  angle  is 
equal  to  the  angle  of  the  incident  ray. 


51 


On  the  other  axis  of  the  mirror,  reflections  follow 


the  rules  of  a  convex  mirror.  That  is,  the  mirror  surface 
has  a  virtual  focus  point  and  all  of  the  reflected  light 
rays  appear  to  come  from  that  focus  point  behind  the  mirror 
surface . 

Reflections  from  other  portions  of  the  cylindrical 
mirror  are  more  complex  than  the  reflections  just 
discussed,  however,  the  overall  effect  is  that  the 
reflection  in  the  cylindrical  mirror  is  compressed  along 
its  main  axis  and  appears  larger  as  the  viewer  moves  closer 
to  the  mirror.  Figure  34  presents  the  view  that  a  camera 
or  viewer  would  see  when  viewing  the  tubing  from  a  position 
out  of  the  plane  of  the  Figures  and  mirror.  This  viewpoint 
is  important  as  it  is  similar  to  the  actual  positioning  of 
the  sensor,  concentrator  and  camera  on  the  solar  thruster. 

Figure  35  presents  the  view  of  the  concentrator  from 
the  viewpoint  of  the  absorber  looking  back  at  the 
concentrator.  The  blue  colored  regions  represent  the 
sunlight  that  might  appear  in  the  concentrator  as  the  focal 
spot  moved  into  correct  focused  alignment.  That  is,  the 
concentrator  should  be  full  of  sunlight  or  almost  full  as 
the  focal  spot  moves  into  the  proper  location  on  the 
absorber . 
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Figure  36  indicates  a  flux  tube  calculation  to  verify 
that  with  perfect  optics,  the  CCD  camera  would  be  able  to 
view  all  of  the  concentrator  by  imaging  the  cylindrical 
mirrors  of  the  absorber.  The  comparison  made  in  Figure  36 
was  between  the  power  received  by  the  camera  by  directly 
viewing  the  Sun  and  by  viewing  the  absorber/concentrator 
system.  Assuming  perfect  optics  it  was  concluded  that  the 
power  received  via  both  paths  was  the  same  so  that  we  were 
guaranteed  that  all  of  the  concentrator  would  be  visible  in 
the  camera  as  well  as  detectable. 

The  new  sensor  system  to  be  described  in  this  chapter 
combines  the  effects  of  the  cylindrical  mirror  and  the  wave 
front  sensor  to  allow  tracking  of  the  solar  focal  spot  on 
the  absorber. 

Sensor  Description 

The  sensor  developed  in  this  section  consists  of  a 
conical  shaped  absorber  made  from  one  guarter  inch  {W') 
metal  tubing.  For  actual  application,  the  tubing  would  be 
made  from  a  high  temperature  metal  such  as  Rhenium  or 
Tungsten.  The  lens  or  optic  system  positions  the  CCD  array 
to  a  virtual  focal  spot  location  behind  the  conical 
absorber.  The  tubes  from  the  conical  absorber  act  like 
cylindrical  mirrors  (toroidal  mirrors)  in  reflecting  the 
sunlight  to  the  CCD.  The  virtual  location  is  important  in 
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applying  the  principles  from  wave  front  sensing,  since  the 
wave  front  sensor  uses  either  apertures  or  lenslets 
projected  onto  the  CCD  array. 

The  lens  system  then  virtually  or  mathematically 
positions  the  array  behind  the  absorber  in  a  manner  similar 
to  the  Shack-Hartmann  system.  The  major  difference  between 
the  method  being  presented  in  this  paper  and  the  Shack- 
Hartmann  system  is  that  the  new  sensor  tries  to  maximize 
the  number  of  wave  fronts  to  determine  focal  spot  location, 
whereas,  the  Shack-Hartmann  system  is  normally  used  to 
optimize  one  particular  wave  front. 

Algorithm 

Each  tube  would  have  its  own  corresponding  subset  of 
virtual  pixels  mapped  directly  on  the  CCD  array.  These 
subsets  would  not  overlap.  Then  each  image  could  be 
processed  by  utilizing  profile  information  along  diagonal 
lines  representing  the  four  quadrants  of  the  circle. 

(Along  the  45  degree  angles,  say)  The  areas  of  maximum 
intensity  could  be  determined  along  each  profile  line.  The 
maximums  should  occur  roughly  where  the  tubes  appear  in  the 
image  as  that  is  where  the  pixel  virtual  images  are 
located.  The  difference  or  gradient  between  these  areas 
should  give  an  indication  of  the  direction  to  the  focal 
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spot.  (Almost  a  centroiding  operation  on  the  maximum  areas 
in  the  image) 

By  knowing  where  the  center  of  the  absorber  is  located 
with  respect  to  the  camera  (a  non-trivial  assumption  as  the 
camera  would  probably  be  mounted  to  one  of  the 
concentrator's  movable  struts),  the  computer  should  be  able 
to  generate  x,  y,  z,  roll,  pitch,  and  yaw  commands  for  the 
hexapod  controller  to  move  the  concentrator  to  a  new 
position  to  provide  better  focus  and  thus  better  heating. 
Figure  37  shows  the  schematic  of  the  proposed  sensor 
solution . 

Schematic  of  Proposed  Solution 


Figure  37  Sensor  and  Camera  Schematic 

For  the  development  of  the  method  presented  in  this 

section,  stainless  steel  H  inch  tubes  were  used  for 
preliminary  imaging,  while  a  copper  conical  absorber  with 
water-cooling  will  be  used  for  refinement  of  the  method. 
Figure  38  shows  the  stainless  steel  tubes  and  Figure  39 
presents  the  water-cooled  conical  shaped  absorber  as  shown 
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in  the  schematic  in  Figure  37. 


The  fluid  to  be  heated  is 


flowed  through  the  absorber  tubing  and  heated  by  the 
concentrated  light. 


Figure  38  STAINLESS  STEEL  TUBES 


Figure  39  COPPER  CONICAL  ABSORBER 


Calculations 

Ray  tracing  using  vectors  was  used  to  determine 
whether  the  image  of  the  7-inch  diameter  concentrator  model 
would  reflect  off  of  the  tubes  and  into  the  aperture  of  the 
CCD  camera.  Figure  40  illustrates  the  situation  and  the 
equation  shown  could  be  used  in  an  iteration  to  find  the 
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angular  extremes  that  would  in  fact  be  imaged  in  the  CCD 
camera.  The  Figure  indicates  a  theta  angle  that  would  be 
double  the  angle  with  respect  to  the  normal  vector  at  the 
surface  of  the  cylinder.  This  is  happening  because  of  the 
way  things  were  drawn  in  Figure  40  and  the  Law  of 
Reflection:  the  angle  of  incidence  is  equal  to  the  angle  of 
reflection . 

A  quick  calculation  using  just  the  rays  from  the 
extremes  of  the  concentrator  model  using  the  Law  of 
Reflection  showed  that  the  CCD  matrix  or  camera  was  able  to 
view  all  of  the  reflection  of  the  concentrator  model  from 
the  tubes. 


Figure  40  VECTOR  CALCULATION  FOR  REFLECTION  ONTO  CCD 

The  field  of  view  of  the  CCD  camera  was  determined  to  be 

1.4  (495  rows)  X  2  (657  columns)  inches  at  a  distance  of  1 

meter  from  the  tubes  and  can. 
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Experiment  3 


A  schematic  of  the  experimental  set  is  shown  in  Figure 
41  and  a  photo  of  the  experimental  equipment  is  shown  in 
Figure  42 . 


Schematic  of  Test  Setup 


CCD  Camera 


Figure  41  TEST  SCHEMATIC 


Figure  42  TEST  SETUP 

The  CCD  camera  used  was  an  ST237  CCD  camera 
manufactured  by  Santa  Barbara  Instrument  Group.  A  657  by 
495  CCD  matrix  with  7.4-micrometer  pixels  characterizes  the 
ST237  CCD  camera.  A  very  small  aperture,  to  gain  depth  of 
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field  and  a  100  mm  lens  form  the  optics  for  the  camera. 

For  the  experiments  discussed  in  this  paper  it  is  focused 
for  objects  located  one  meter  away  from  the  lens. 

The  concentrator  model,  see  Figures  9X4,  10X4,  and 
11X4,  consisted  of  a  7-inch  diagonal  circle  of  white 
surrounded  by  a  black  mask.  At  the  distances  shown  in 
Figure  7X4,  the  concentrator  half-angle  model  was 
approximately  7  inches  in  diameter.  A  standard  overhead 
projector  was  used  to  supply  light  to  the  concentrator 
model  using  a  back  projection  method. 

Images  of  the  stainless  steel  tubes  were  taken  to 
discern  whether  an  image  of  the  full  view  concentrator 
model  could  be  viewed  in  the  tube.  Several  different 
exposure  times  were  used  as  well  as  different  concentrator 
patterns  from  full  white  circles  to  circles  with  varying 
amounts  of  area  blackened  to  represent  various  stages  of 
misalignment . 

The  experiment  proceeded  by  first  using  one  vertical 
tube,  and  then  two  vertical  tubes.  The  three  patterns  used 
are  shown  in  Figures  43  through  45.  Figure  43  shows  the 
full  concentrator  pattern.  Figure  44  shows  pattern  one 
representing  a  concentrator  configuration  that  would  be 
considered  not  on  focus,  and  Figure  45  shows  pattern  number 


59 


two  with  even  less  white  in  the  concentrator  area 
representing  less  focus  than  pattern  one. 


Figure  43  9X4  FULL  CONCENTRATOR  MODEL 


Figure  44  PATTERN  ONE  MODEL 


Figure  45  PATTERN  TWO  MODEL 

ONE  TUBE 

The  single  tube  images  indicated  that  our  camera  might 
not  have  enough  resolution  to  be  able  to  discern  between 
concentrator  pattern  images.  Figure  46  shows  the  image  of 
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one  tube  with  the  full  concentrator  pattern,  and  Figure  47 
shows  the  profile  through  the  image. 

The  images  and  profiles  of  the  images  indicated  very 
little  differences  between  the  patterns.  The  single  tube 
experiment  was  supposed  to  show  a  complete,  although 
distorted  image  of  the  concentrator  light.  At  the 
resolution  of  the  camera,  7 . 4um  pixels  focused  at  1  meter 
away  from  the  cylinder  did  not  allow  a  complete  picture  of 
the  concentrator  model.  This  experiment  may  need  to  be 
repeated  after  a  camera  with  higher  resolution  is  acquired. 


Figure  4  6  ONE  TUBE  FULL  CONCENTRATOR 


Distance  along  profile 


Figure  47  PROFILE  FOR  FIGURE  46 

TWO  TUBES 
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The  images  from  the  two-tube  experiment  are  shown  in 
Figures  48-50  and  their  profiles  are  shown  in  Figures  51- 
53.  When  each  pattern's  image  and  profile  are  examined, 
the  profiles  showed  a  distinct  lowering  of  the  peak 
intensities  seen  at  each  tube  location.  Pattern  number  2 
showed  the  smallest  peaks  compared  to  the  full  concentrator 
model  and  to  pattern  number  1.  Since  pattern  2  had  more 
black  area  in  the  white  concentrator  area,  it  was  expected 
to  produce  less  intensity  on  the  tubes.  These  results 
indicated  that  the  concept  of  the  new  sensor  would  work  as 
thought.  However,  since  the  patterns  were  not  directly 
related  to  actual  misalignment  of  the  concentrator,  no 
actual  position  information  could  be  derived  from  the 
profile  information. 


Figure  48  Two  Tubes,  Full  Concentrator 
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Figure  49  Two  Tubes  Pattern  1 


Figure  50  Two  Tubes,  Pattern  2 


Figure  51  Profile  Two  Tubes  Full  Concentrator 
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Figure  52  Profile  Two  Tube  Pattern  1 
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Figure  53  Profile  Two  Tube  Pattern  2 

Summary  Experiment  3 

From  the  experiments  conducted  for  this  section,  the 
concept  of  using  the  cylindrical  mirrors  of  the  tubing  for 
the  absorber  as  a  wave  front  sensor  appears  to  be  feasible 
and  correct.  The  method  of  using  one  tube  to  look  at  the 
concentrator  model  revealed  the  fact  that  the  CCD  camera 
being  used  does  not  have  enough  resolution  to  image  the 
concentrator  reflection  at  1  meter.  The  two-tube 
experiment  showed  a  very  promising  pattern  in  that  it 
showed  a  decrease  in  intensity  in  conjunction  with  a 
reduction  in  the  white  areas  Ought  from  the  full 
concentrator  QConcentrator  to  the  other  two  patterns.  Since 
the  patterns  did  not  represent  actual  misalignment,  further 
work  would  be  needed  to  make  suitable  patterns  to  represent 
misalignment . 

Subsequent  work  could  focus  on  improving  the 
concentrator  models  by  making  them  conform  to  actual 
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misalignment  and  to  making  the  black  areas  able  to  block 
more  light.  A  possible  error  source  in  the  images  may  be 
IR  light  from  the  overhead  projector  transmitted  through 
the  black  paper  and  causing  the  intensities  to  depend  more 
on  the  IR  instead  of  just  the  white  light.  Also,  since  the 
projector  was  not  completely  black  around  the  areas 
backlit,  scattered  light  may  have  entered  the  CCD  and 
caused  intensity  changes. 

A  quick  way  to  make  the  black  outlines  darker  and 
opague  to  IR  would  be  to  either  back  the  model  with 
cardboard  or  make  the  patterns  out  of  Aluminum  foil.  A 
cardboard  mask  that  masked  all  but  the  circular  area  of  the 
concentrator  could  be  used  to  improve  the  extraneous  light 
reaching  the  CCD  camera. 

4.4  Concentrator  Model 

The  concentrator  model,  to  model  actual  misalignment, 
would  have  to  be  constructed  such  that  the  dark  areas 
appear  at  the  edges  of  the  concentrator  first  as  the 
concentrator  moves  out  of  focus.  At  that  point  the 
concentrator  would  have  a  white  center  with  black  edges. 
Modeling  and  simulation  of  the  concentrator' s  ability  to 
capture  sunlight  was  used  to  determine  angles  and  distances 
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needed  in  the  imaging  of  the  absorber  sensor  using  an 
actual  concentrator. 

The  copper  coiled  absorber  sensor  would  become  the 
primary  target  for  images  of  the  concentrator  model.  An 
actual  1X2  meter  concentrator  was  available  to  perform 
more  analysis  of  the  image  formed  in  the  concentrator  when 
off  focus.  Finally,  a  camera  with  more  resolution  was 
obtained  to  allow  more  experimentation  with  the  tubes  to 
investigate  the  full  image  of  the  concentrator  providing 
more  insight  to  the  functioning  of  the  cylindrical  mirrors 
and  wave  front  sensor. 

A  JAVA  program,  Concentratorv3 . j ava,  was  written  to 
model  the  sunlight  collected  by  a  paraboloid  concentrator. 
The  program  listing  is  in  Appendix  B.  The  concentrator 
model  program  traced  light  rays  from  the  target  sensor  (the 
absorber)  to  the  Sun.  The  diagram  in  Figure  54  shows  the 
vector  construction  used  in  the  programming  of 
Concentratorv3 . j  ava . 
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Figure  54  Vector  Construction  For  JAVA  Model 

The  basic  operation  of  the  program  is  that  a  square 
grid  is  built  at  the  eye  position  which  then  views  the 
concentrator.  The  concentrator  view  is  formed  on  the  grid 
as  a  projection  of  the  concentrator  surface  onto  the  grid 
plane  pointed  to  by  eoo  and  positioned  in  space  using  angle 
from  z  axis  and  c,  cone  half  angle  using  the  coordinate 
system  shown  in  Figure  55.  The  position  of  the  "eye"  or 
viewpoint  vector,  initially  located  at  the  focus  of  the 

paraboloid,  is  located  in  space  by  the  vector^  .  The 
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VectorX,  locates  each  of  the  incremental  mirrors  or  pixels 
that  make  up  the  surface  of  the  concentrator. 

As  the  program  is  run,  ei,j  is  incremented  through  in 
each  direction  filling  in  the  grid  of  the  eye  or  CCD  for 
each  vector  ei,j.  A  grid  point  on  the  CCD  is  filled  if  the 
angle  between  the  reflected  ray,  vector  T  in  Figure  54, 
and  vector  R  is  less  than  the  angle  of  the  Sun,  0.5 
degrees.  Each  run  of  the  program  collects  information 
regarding  one  location  on  the  target  or  absorber  plane. 


Vector  From  CCD  or 
Targetto  center  of 

Angle  of  if  From  the 
zAxis 


Vector  From  CCD  or 

Concentrator  under 
Consideration 


Vectorfrom  origin 
to  CCD  location 


Image  Coordinates 
u  =  Au*i*e* 
v  =  Av*j*e2 


Figure  55  Coordinate  System  For  Target  Plane 

The  datum  returned  from  the  program  was  1  if  the  ray 

could  be  traced  back  to  the  Sun  and  a  0  if  the  ray  could 
not  be  traced  back  to  the  Sun.  Three  Figures  are  included 
to  show  the  effects  of  being  off  focus  by  very  small 
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offsets  to  S  and  to  various  changes  to  and  z.  The 
Figures  are  Figures  56,  57,  and  58. 


Thetaz  =  0.5  degrees,  Phiz  =  0  degrees,  Eye  0.01  m 


Figure  56  Thetaz  =  0.5  Degrees,  Phiz  =  0,  Eye  =  0.01m 


Thetaz  =  0.5  degrees  Phiz  =  10  degrees  Eye  =  0.01  m 


Figure  57  Thetaz  =  0.5  Degrees,  Phiz  =  10  Degrees,  Eye  =  0.01m 
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Thetaz  =  0.5  degrees.  Phiz  =  20  degrees,  Eye  0.01  m 


Figure  58  Thetaz  =  0.5  Degrees,  Phiz  =  20  Degrees,  Eye  =  0.01m 

To  produce  a  full  picture  of  the  light  in  the  target 
area,  one  would  have  to  either  run  the  program  multiple 
times  and  save  the  information  for  collating  or  modify  the 
program  to  run  multiple  view  or  eye  points  saving  the  final 
data  for  plotting.  The  way  the  program  currently  functions 
was  sufficient  to  determine  that  the  control  environment 
for  the  focal  spot  location  was  within  one  centimeter  of 
the  designed  focal  spot.  That  is,  very  little  perturbation 
of  the  pointing  of  the  focal  spot  reduced  the  area  of  the 
sunlight  being  collected  and  focused  on  the  target  or 
absorber  to  the  point  that  no  light  was  on  the  absorber 
sensor . 
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4.5  Algorithm  Development 


Control  of  the  position  of  the  solar  concentrators  to 
ensure  proper  heating  of  the  thruster  is  paramount  in  STP. 
Prior  work  settled  on  an  arbitrary  control  requirement  of 
0.1  degree  in  angular  position  and  0.1  inch  in 
translation5™ . 

Continuing  from  the  definition  of  Isp  from  the  section 
on  solar  thermal  propulsion  background,  Isp  is  the  thrust 
divided  by  the  propellant  expelled  in  the  exhaust. 


F 


1 


I  sp 

m  m 


dp 

dt 


=  —mve  =ve 
m 


(4: 


In  4,  ve  is  the  exhaust  velocity.  From  4,  the  expression 
for  the  power  developed  by  the  exhaust  stream  of  the 
thruster  is  shown  in  5. 


P  =  F  mv  =  —  ml  2 
2  p 


(5; 


Equation  5  can  be  related  to  concentrator  power  input  Pjn  to 
the  thruster  and  the  thruster  loss  mechanism  r\t  as  follows 
in  equation  6. 

P  =V,Pin  (6) 
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For  the  concentrator  portion  of  Pjn  ,  the  equation  for 


W  W 

the  power  is  in  equation  7  with  I0  =  1350 — -±50 — -  at  Earth 

in  “  m 


with  aperture  area  Ac  . 


P  =  hA  (7) 

The  power  transferred  to  the  thruster  is  related  to 
the  intensity  of  the  light  collected  by  the  concentrators. 
The  power  in  equation  7  is  modified  as  the  thruster 
reflects  some  power  away  from  the  thruster  and  radiates 
energy  at  the  thruster's  black  body  radiation  temperature. 
The  power  into  the  thruster  is  then  shown  in  equation  8. 

(8) 

Where  £  is  the  emissivity  of  absorber  area  A,  cris  the 
Maxwell-  Boltzmann  constant,  Iin  is  the  input  power  density 
from  concentrator  to  absorber  A,  and  T  is  the  temperature 
of  the  absorber  and  propellant.  Equation  8  is  a  simple  one 
dimensional  expression  for  the  power  into  the  thruster. 

The  assumptions  for  equation  8  are  that  the  thruster 
has  a  uniform  temperature  over  A,  the  temperature  over  A  is 
the  same  for  the  propellant,  and  £  is  equal  to  1  minus  the 
reflectivity.  Otherwise  equation  8  becomes  an  integral 
over  the  flow  path.  In  equation  8,  Iin  is  found  by  relating 
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the  solar  power  density  with  effects  of  the  optical  path 
through  the  concentrator  as  shown  in  9. 

I  in  =*CCrI0  (9) 

Where  tc  is  the  transmissivity  of  the  concentrator,  Cr  =  -j 
and  Ac  is  the  area  of  the  concentrator  aperture  and  A  is 
the  area  of  the  thruster/absorber ,  and  I0  is  the  solar  power 
density  at  Earth  orbit.  From  equations  9  and  6,  the  power 
at  the  exhaust  stream  is  shown  in  equation  10. 

pe  =rJt£A^cCrh  -ctr4J  (10) 

Equations  4  and  10  may  be  used  to  find  thrust  F  as  a 
function  of  I  ,  T,  and  I0  . 

Relating  the  power  to  the  exhaust  power,  Pe  to  the 
misalignment  angle  of  the  concentrator  is  a  problem  solved 
numerically  using  the  concentrator  code  written  by  Dr. 

Holmes  called  offaxis.  Figure  59  shows  the  result  for 
several  angles  of  misalignment.  Using  values  from  Figure 
59,  Figure  60  plots  the  intensity  versus  misalignment 
angle.  Figure  60  shows  that  the  power  density  decreases 
considerably  at  0.2  degrees  of  misalignment. 
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Figure  59  Peak  Power  Density  vs .  Radius  Off  Focus 

The  peak  power  at  0.2  degrees  of  misalignment  is 
approximately  2  dB  down  and  indicates  a  likely  cutoff  point 
for  this  discussion.  Therefore,  most  of  the  work  presented 
in  this  article  pursued  a  0.1  degree  of  misalignment  target 
value  for  maximum  power  transfer  to  the  thruster. 
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Intensity  Vs  Solar  Misalignment 


Figure  60  Intensity  vs .  Solar  Misalignment 

Section  4.3  presented  the  idea  of  using  the  absorber 

of  the  thruster  as  a  sensor  for  tracking  of  the  focal  spot 
on  the  absorber.  This  section  presents  two  algorithms  used 
in  determining  the  location  of  the  focal  spot  using  the 
absorber  as  the  primary  sensor.  Both  algorithms  utilize 
concepts  developed  from  studies  of  the  Shack-Hartmann  wave 
front  sensor.  The  first  algorithm  described  uses 
correlation  and  masking  for  locating  and  tracking  the  focal 
spot.  The  second  algorithm  simplifies  the  method  in 
algorithm  one  to  using  area  moments  for  tracking  sunlight. 
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System  Description: 


Since  section  4.2  presents  the  hardware  being  used  for 
focal  spot  tracking  in  detail,  a  brief  summary  of  the 
components  will  now  be  presented  for  use  in  this  section. 
The  schematic  of  the  hardware  is  shown  in  Figure  71  and 
consists  of  a  CCD  camera  and  a  conical  shaped  tubular 
absorber . 

The  new  sensor  and  algorithms  are  based  on  the  Shack- 
Hartmann  wave  front  sensor  with  the  holes  and  lenslets 
replaced  by  the  cylindrical  mirrors  of  the  current  system. 
Sunlight  reflecting  on  the  absorber  tubing  is  imaged  by  the 
camera  for  analysis  instead  of  imaging  the  array  of  holes 
or  lenslets  of  the  Shack-Hartmann  device. 

The  thruster  tubes  serve  as  high  temperature  mirrors 
or  "lenslets."  The  absorber  takes  the  heat  from  the 
sunlight  and  conducts  it  to  the  Hydrogen  gas  for  use  in 
providing  thrust.  Differences  in  location  of  light  and  how 
much  light  is  reflected  from  the  tubes  are  used  to 
determine  the  focal  spot  location  error  for  use  in  a 
feedback  control  system. 
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Two  or  possibly  three  different  algorithms  will  be 
used  to  locate  the  focal  spot  error  from  the  images  of  the 
tubes.  Different  algorithms  were  deemed  necessary  because 
of  the  0.1  degree  arbitrary  control  requirement  and  not 
having  one  universal  algorithm  that  covered  the  three 
regimes  of  focus  (sunlight  more  than  5  diameters  off 
focus) ,  intermediate  focus  (approximately  5-3  diameters  off 
focus) ,  and  fine  focus  (2-1  diameters  off  focus) . 

Algorithms: 

METHOD  ONE: 

The  image  processing  portion  of  the  focus  control 
system  was  developed  along  two  paths.  The  first  method  to 
be  discussed  is  the  more  general  method  using  correlation, 
while  the  second  method  is  less  general  based  on  area 
moments  and  area  centroids  similar  to  one  version  of  the 
Shack-Hartmann  analysis.  The  flow  chart  in  Figure  61  shows 
the  sequence  to  be  implemented  for  solar  concentrator  focus 
control . 
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Figure  61  Flow  Chart  of  Sequence 


Figure  62  illustrates  the  diagram  of  the  system  after 
deployment  and  before  coarse  focus.  The  method  for  coarse 
focus,  which  is  not  discussed  in  this  dissertation,  works 

on  minimizing  §  to  get  to  the  situation  shown  in  Figure  63 
The  situation  is  then  the  initial  position  for  the  fine 
focal  controller,  the  focal  spot  is  located  in  the  general 
area  of  the  sensor  but  not  on  the  optimized  spot  and  the 
system  is  not  thrusting.  Thus,  after  deployment,  the 
camera  is  pointed  at  the  absorber  and  the  concentrator  is 
positioned  facing  the  Sun.  This  situation  is  then  defined 
as  the  initial  position  for  the  coarse  control  algorithm 
shown  in  Figure  62 . 
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With  the  addition  of  a  stored  image  of  the  absorber  in 


the  on-Sun  and  on-focus  designed  and  geometric 
configuration,  the  control  algorithms  begin. 


Figure  62  Starting  Coordinates  For  Coarse  Focus 


Figure  63  Fine  Focus  Start  Coordinates 

At  the  start  of  the  control  program,  after  coarse 
focus,  the  control  computer  has  the  initial  geometric 
designed  configuration  and  initial  image  stored  in  memory. 
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These  initial  conditions  also  define  an  initial  coordinate 


system  as  displayed  in  Figure  63.  Also  stored  in  the 
control  computer  are  masks  to  be  used  in  computing  where 
the  focal  spot  is  currently  located  each  iteration. 

The  correlation  or  the  lack  of  correlation  at  the 
masked  areas  indicates  the  current  location  of  the  spot  and 
indicates  the  direction  of  travel  required  to  reduce  the 
mismatch.  The  overall  process  for  this  control  algorithm 
follows  image  correlation  or  scene  matching  from  computer 
pattern  recognition*™1,  XV11,  xvi;l;l  .  Both  procedures  follow  a 
similar  mathematical  development  to  that  of  correlation 
between  masked  areas  of  the  test  image  with  the  image  of 
the  on-focus  image.  Optimization  would  be  utilized  to  find 
the  optimal  and  correct  focal  spot  location. 

From  Figure  63 ,  rerror  ,  p  ,  X(d),  and  nerror  are  the  primary 
vectors  in  the  development  of  the  algorithms.  When  the 
system  is  deployed  on  orbit,  the  initial  position  of  the 
concentrator,  camera,  and  thruster  determine  the  initial 
origin  and  coordinate  system.  ferror  is  the  vector  that 
determines  the  location  of  the  image  in  the  original 
coordinate  system  for  the  focus  spot.  nerror  is  the  vector 
normal  to  the  plane  at  the  current  image  location  atrerror  . 
p  ,  is  the  location  of  the  camera  or  "eye-position"  with 
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respect  to  the  original  coordinate  system.  Finally,  the 
vector  x(g)  represents  areas  of  the  images  that  are  regions 
of  interest  for  locating  the  focal  spot.  The  regions  of 
interest  in  the  case  of  these  algorithms  would  be  the 
thruster  itself  or  the  individual  rings  of  the  absorber. 

The  required  accuracy  of  focal  spot  location  determines 
whether  the  thruster  as  a  whole  would  be  the  region  of 
interest  or  the  individual  rings. 

For  the  imaging  portion  of  the  process,  the  vector  x(g) 
in  the  proceeding  paragraph  would  have  to  be  projected  onto 
the  CCD  camera  image  plane.  The  method  used  for  this 
operation  is  based  on  pinhole  optics  and  perspective 
projection.  The  CCD  camera  utilizes  a  pinhole-like  or  small 
aperture  for  its  almost  infinite  depth  of  field.  A  lens  is 
positioned  just  behind  the  aperture  to  focus  the  aperture 
image  onto  the  CCD  matrix  in  the  camera.  With  these 

optical  parameters  the  projection  of  x(g)  onto  the  CCD 
matrix  with  coordinates  (a,b)  is  determined  using  the 
following  equations. 

a{6)  =  a  *  (jx(0)  +  r error  -  p }  C  )(x(0)  +  rerror  -  p }  e3  (11) 

b{6 )  =  /3  *  ((x(0)  +  Terror  ~p)e2  )^(0)  +  V  error  P  ^  £3  (12) 
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a  and/?  in  equations  11  and  12  are  scaling  factors 
determined  by  the  CCD  pixel  size.  The  e  vectors  are  unit 
vectors  located  at  p  defining  a  coordinate  system  at  p  for 
the  camera. 

Now  that  the  image  and  system  coordinates  and  vectors 
have  been  defined,  the  correlation  function  is  derived. 

The  image  taken  by  the  CCD  will  be  denoted  I (a,  b)  and  is  a 
function  of  ( a  ,  /?  )  as  expected.  W  ( a  ,  /?  )  will  represent 
the  image  or  sub-image  that  we  are  trying  to  match.  For 
example,  W  (a,  /?  )  would  represent  the  ideal  on-focus  image 

or  it  could  represent  only  a  portion  of  the  ideal  image 
such  as  specific  rings  of  the  sensor.  The  L2  norm  will  then 
be  utilized  to  calculate  the  minimum  distance  between  W  (a, 
/?  )  and  the  current  image . 

E(m,  n)  =  ^  Jj=i  [/ O'  +  m  - 1,  j  +  n  - 1)  -  W(1, 7')}  )  (13) 

Where  E(m,n)  is  the  measure  of  distance  at  coordinates 

(m,  n) 

Equation  13  also  assumes  that  W  is  smaller  than  the 
image  I,  but  that  does  not  preclude  using  a  W  that  is  equal 
to  I  in  size.  Squaring  and  expanding  both  sides  of 
equation  13  while  using  the  assumption  that  the  energies  in 
W  and  I  are  constant  and  very  small,  the  correlation  of  W 
and  E  results  as  shown  in  equation  14. 
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(14) 


E2  (m,  ri)  =  -2*  /(/  +  m  - 1,  j  +  n  - 1)*  W (i,  J ) 

Multiplying  equation  14  by  -1/2  generates  the  correlation 
R(m,n)  at  (m,n)  as  shown  in  equation  15. 

R(m,  n )  =  ^  ^  I(i  +  m  -  1,  j  +  n  - 1)  *  W (/,  j)  (15) 

Finding  the  maximum  value  for  R(m,n)  will  determine 
the  best  location  in  the  image  that  matches  W.  Finding  the 
X,Y  location  of  the  maximum  value  for  R(m,n)  and 
determining  the  difference  from  that  location  to  the 
location  of  the  autocorrelation  of  the  on-focus  image 
determines  which  direction  to  move  the  image  to  match  the 
maximum  correlation.  Once  the  difference  between  the 
autocorrelation  of  the  on-focus  image  location  and  the 
point  of  maximum  correlation  is  found,  the  control  system 
would  move  the  concentrators  to  reduce  this  difference. 
Equation  15,  correlation  was  thought  to  be  sufficient  for 
determining  the  location  of  the  focus. 

MASKS : 

Several  different  masks  were  generated,  using  a  free 
hand  drawing  method,  to  help  determine  which  masks  were 
best.  Best  in  this  usage  means  that  the  masks  were  best 
suited  to  determining  the  offset  distance  from  the  focus 
location.  Several  masks  that  were  used  with  the  tubular 
absorber  are  shown  in  Figures  64,  65,  and  66.  Similar  masks 
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were  used  with  the  loop  of  tubing  that  replaced  the 
absorber  for  testing. 

The  difference  in  the  masks  was  which  ring  of  the 
absorber  was  used  in  preparing  the  mask.  For  the  outer 
ring  mask,  the  outermost  ring  of  the  absorber  was  used  as 
the  masking  area  and  thus  enclosed  a  bigger  area  of  the 
image  than  that  covered  by  the  inner  ring  mask.  So 
although  the  masks  in  Figures  64  and  65  look  alike,  they, 
in  fact,  cover  very  different  regions  of  the  images. 

The  mask  in  Figure  66  was  thought  to  be  able  to 
isolate  a  portion  of  the  outer  ring  for  consideration.  In 
actual  practice,  it  matched  too  many  intensity  levels  to  be 
of  use,  especially  once  the  light  got  less  than  4  diameters 
off  of  focus.  In  other  words.  Figure  66  could  not  be  used 
to  discriminate  below  about  4  diameters  as  its  correlation 
was  smeared  and  did  not  produce  impulse-like  behavior  at 
correlation . 


Figure  64  Outer  Ring  Mask 
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Figure  65  Inner  Ring  Mask 


Figure  66  "Banana"  Shaped  Mask 


METHOD  TWO: 

The  second  method  for  determining  focal  spot  location 
is  by  measuring  light  area  moments  in  the  image.  This 
method  is  derived  based  on  a  modification  of  the  Shack- 
Hartmann  wave  front  sensor  using  cylindrical  mirrors 
instead  of  lenslets  to  register  light  areasxlx.  Images  of 
the  concentrator  can  be  seen  in  each  coil  of  the  absorber. 
Whether  a  particular  image  on  the  tube,  denoted  as  a  "tube 
blob"  is  bright  or  dim  depends  on  how  close  the  tube  is  to 
the  focal  point.  Image  area  moments  are  then  used  to 
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quantify  focus  for  individual  "tube  blobs." 


The  derivation 


also  follows  from  a  discussion  of  area  moments  and 
centroids  from  mechanics™,™1.  An  image  is  acquired  via 
the  CCD  camera  and  a  threshold  is  applied  to  the  image  for 
further  processing. 


T ube  "Blobs”  can  appear  anywhere 
along  or  on  the  tube. 


Plate  “Blobs”  appear  on  the  plate 
in  the  center  of  the  tube. 


Figure  67  Generic  Light  Patterns  in  Image  of  Sensor 


Figure  67  shows  a  generic  schematic  of  light  blobs 
from  an  image  of  the  sensor.  As  explained  earlier,  each 
tube  "blob"  is  an  actual  image  of  the  light  from  the 
concentrator.  The  shape  and  intensity  area  of  each  tube 
blob  is  related  to  the  amount  of  light  in  the  concentrator 
being  concentrated  to  that  area  of  the  tube.  Calculating 
the  area  centroid  of  the  images  in  the  tubes  determines 
where  the  focal  spot  is  located  on  this  iteration  of  the 
algorithm  and  knowing  where  the  centroid  of  the  focal  spot 
is  located  for  this  iteration  determines  which  direction 
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the  control  algorithm  needs  to  move  the  concentrator  to 


move  the  centroid  to  the  center  of  the  sensor. 


Figure  68  Light  in  Image  From  CCD 

Figure  68  shows  the  situation  that  might  occur  once 
the  CCD  acquires  an  image.  A  coordinate  system  for  the 
image  is  setup  as  shown  in  Figure  68.  The  image  is  then 


processed  by  first 

thresholding 

the 

image  and 

then 

calculating  the  area 

moments  for 

each 

light  area. 

The 

threshold  is  set,  at 

least  in  the 

ground  experiments 

,  by 

visually  making  sure  that  the  lighted  areas  in  the  sensor 
tubes  do  not  overlap  into  the  area  of  the  secondary 
concentrator  and  contain  approximately  95%  of  available 
light  in  the  images  reflected  from  the  concentrator  when  on 
focus.  On  flight  another  method  of  preventing  overlap 
would  be  needed,  possibly  a  measured  area  of  the  absorber 
beyond  which  the  controller  would  ignore.  Once  the  image 
has  been  thresholded,  the  area  moments  are  calculated  and 
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an  average  or  overall  centroid  is  calculated  as  the 


location  of  the  focal  spot  on  the  sensor. 


Figure  69  Calculation  Schematic 

Figure  69  is  an  ideal  schematic  of  the  calculation  for 
one  of  the  ellipse  shaped  light  areas.  The  light  area  is 
located  at  r2  and  the  area  calculated  as  (*a*b)/4,  where 
(a)  is  the  major  axis  of  the  ellipse  and  (b)  is  the  minor 
axis  of  the  ellipse.  Position  vector  r2  is  the  vector  from 
the  origin  to  the  center  of  the  ellipse. 
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Figure  70  Centroid  Calculation 

Finally,  Figure  70  shows  the  calculations  shown  in 
Figures  68  and  69  combined  to  calculate  the  overall  area 

centroid  as  noted  by  {X,Y).  This  position  is  the 
calculated  focal  spot  location.  As  shown  in  Figure  70,  the 
image  of  the  central  flat  plate  will  have  a  significant 
"blob"  of  light  located  at  the  center  of  the  sensor  as  the 
control  described  in  this  dissertation  assumes  a  coarse 
focus  is  obtained  before  the  methods  of  this  paper  are  used 


for 

fine  focus. 

If 

the 

"blob"  of 

light 

in  the  center 

of 

the 

image 

is 

too 

big 

for  the 

area 

centroids  to 

be 

calculated. 

it 

can 

be 

suppressed 

in  the  calculation 

of 

centroids  to  obtain  locations  that  are  very  close  to  focus. 

Experiment  4 

Figure  71  depicts  the  setup  used  for  collecting  images 
for  analysis.  The  source  of  light  used  was  a  3  inch 
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diameter  red  LED  taillight.  The  size  of  the  light  gives  an 
extended  source  a  little  over  half  of  the  Sun'’  s  diameter  as 
viewed  from  the  concentrator.  For  the  setup  shown  in 
Figure  71,  the  size  of  light  that  corresponds  to  Earth's 
orbit  is  approximately  5  inches  in  diameter.  SRS 
Technology' s  1X2  meter  elliptical  concentrator  was  used  to 
focus  light  on  the  thruster/absorber  and  secondary 
concentrator  as  well  as  the  later  imaging  of  the  single 
loop  of  tubing  utilized  in  the  study  of  misalignments 
closer  than  2-3  diameters  off  of  focus. 

Images  were  obtained  using  a  STX-10XE  CCD  camera  from 
Santa  Barbara  Instrument  Group  (SBIG) .  This  camera 
provides  a  2184  X  1472  pixel  6.8  CCD  for  imaging.  CCDOPS 
from  SBIG  was  used  on  a  notebook  computer  to  acguire  images 
at  0.5  sec  exposure  times.  The  Sun  simulator  was  connected 
to  a  tripod  which  allowed  the  Sun  simulator  to  move  around 
at  the  second  focal  spot  of  the  concentrator  providing  the 
misalignment  of  the  Sun  and  the  concentrator  for  study. 
Images  were  then  taken  with  the  Sun  simulator  at  various 
locations  from  focus  in  1  diameter  increments.  The  maximum 
distances  that  were  obtainable  in  the  physical  layout  shown 
in  Figure  71  was  3-4  diameters  off  of  focus  in  the  X 
direction  and  2-3  diameters  in  the  Z  direction.  Movement 
in  the  y  direction  was  not  checked  in  these  experiments. 
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Figure  72  and  Figure  73  show  the  actual  equipment  used  in 


the  experiments.  The  assigned  coordinate  axes  for  the 
light  source  are  shown  in  Figure  74  and  in  Figure  75  for 
the  thruster  in  these  experiments.  These  coordinate 
systems  become  more  important  for  control  of  the  hexapod 
used  for  moving  the  concentrator. 


Test  Apparatus 


Figure  71  Test  Setup 


Figure  72  Concentrator  From  Source  End 
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Figure  73  Divergent  Source 


Y 


X 


Figure  74  Top  View  Coordinate  System  for  Light  Source 


Y 


X 


Figure  75  Top  View  Coordinates  at  Thruster 

ImageJ  and  Image  Analysis™" 

The  ImageJ  image  analysis  and  processing  package  from 
the  National  Institutes  of  Health  was  used  to  process  the 
images  taken  in  the  experiments.  The  program  provides  a 
Graphical  User  Interface  (GUI)  and  extensive  image  analysis 
commands.  ImageJ  is  also  extendable  by  using  either  macros 
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or  ImageJ  plugins.  Plugins  are  basically  Java  classes  that 
implement  one  of  several  interfaces  in  ImageJ  or  they  are 
Java  classes  performing  specific  operations  within  ImageJ. 

Method  1 

ImageJ  performs  correlation  on  images  using  the 
conjugate  multiplication  of  the  Fourier  transforms  of  two 
images  to  be  correlated.  Masks  were  utilized  to  reduce  the 
size  of  images  to  include  only  the  absorber  area  as  the 
region  of  interest.  Masking  also  helped  the  processing 
time  required  for  the  calculation  of  the  FFTs  as  the  images 
were  reduced  from  2184  X  1472  down  to  a  more  manageable  512 
X  512  while  still  maintaining  the  correct  region  of 
interest.  Results  from  the  processing  of  various  offsets 
from  on-focus  are  presented  in  Figures  76  through  79. 

Figure  76  shows  the  autocorrelation  result  for  the  on-focus 
image.  It  has  a  very  clear  "dot"  or  impulse  at  the  center 
of  the  image.  This  result  was  expected  as  the 
autocorrelation  of  an  image  with  itself  would  be  an  impulse 
at  offset  (0,0) . 

Figure  77  shows  the  result  for  the  cross-correlation 
between  the  images  for  up  1  diameter  from  on-focus  of  the 
solar  simulator  and  the  on-focus.  In  this  instance,  the 
correlation  is  smeared  and  basically  located  at  the  center 
of  the  image  indicating  that  the  1  diameter  misalignment 
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has  not  been  resolved  from  the  on-focus  situation  in  which 


the  maximum  value  is  located  in  the  center  of  the  image. 

Figure  78  shows  a  much  better  situation  for  the 
correlation  technique.  This  Figure  shows  the  correlation 
between  a  4  diameter  image  and  the  on-focus  image.  This 
Figure  also  shows  the  result  for  a  mask  similar  to  the  one 
shown  in  Figure  66.  The  mask  used  for  Figure  78  has  other 
problems  at  a  smaller  offset  as  the  numbers  of  correlations 
that  may  occur  in  the  image  get  larger  as  the  focus  spot 
gets  closer  to  the  on-focus  location. 


Figure  76  Autocorrelation 
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Figure  78  Correlation  For  a  4  Diameter  Offset 


Figure  79  Correlation  For  a  2  Diameter  Offset 


The  correlation  shown  in  Figure  79  illustrates  the 
problems  involved  when  using  the  correlation  function  at 
offsets  below  about  3  diameters.  The  correlation  in  Figure 
79  shows  that  an  offset  exists;  however,  the  offset  itself 
cannot  be  determined  from  this  image.  The  cutoff  for  using 
masks  like  those  in  Figures  64  and  65  using  correlation 
defined  by  Equation  15  appears  to  be  between  2  and  3 
diameters  of  offset. 

However,  at  this  point  it  was  noted  that  using  only 
Equation  15  lead  to  problems  as  it  would  find  any  and  all 
correlations  based  on  the  maximum  intensities  present,  with 
little  regard  for  the  variations  in  intensities  present  in 


95 


an  area  of  the  image,  instead  of  the  correct  correlation 
value.  That  is.  Equation  15  would  pick  an  area  in  I(m,n) 
that  in  general  had  very  high  intensities  instead  of 
actually  the  point  of  highest  correlation.  That  problem  is 
eliminated  by  dividing  Equation  15  by  the  total  energy  in 


I  (m,  n)  , 


1 


l  j  i  \ 

^  I1  (i  +  m  - 1  ,j  +  n  -  l) 


Furthermore,  Equation  15 


can  be  normalized  to  the  range  (-1,1)  by  dividing  by  the 
energy  in  W(m,n) .  It  was  not  necessary  to  divide  by  the 
energy  of  W(m,n)  for  the  case  being  discussed  in  this 
dissertation  as  it  was  not  important  to  get  values  from  (- 
1,  1) .  Another  form  for  Equation  15  is  known  as  the  phase 

only  correlation  that  basically  divides  W(m,n)  by  the 
magnitude  of  W(m,n),  which  normalizes  the  FFT  of  W 
providing  only  the  phase  information  which  provides  a  more 
impulse-like  function  for  the  correlation.  As  a  quick 
approximation  of  the  phase  only  correlation,  a  high  pass 
filtering  was  tried  since  the  phase  only  equation  is 
similar  to  Equation  15  but  with  the  addition  of  dividing  by 
the  magnitude  of  W.  Dividing  by  the  magnitude  of  W,  with  W 
having  low  pass  frequency  patterns,  is  approximately  high 
pass.  However,  phase  only  correlation  was  abandoned 
because  it  did  not  handle  rotations  of  the  region  of 
interest  satisfactorily. 
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Method  2 

ImageJ  was  also  utilized  for  the  analysis  outlined  in 
method  two  of  this  paper.  The  particular  command  used  was 
the  "Analyze  Particles"  of  the  Analyze  menu.  To  prepare 
the  images  for  analysis,  each  image  had  to  have  a  threshold 
applied  before  "Analyze  Particles"  would  work.  For  the 
experiments  in  this  paper  the  automatic  threshold  was  used 
for  all  images.  It  was  felt  that  auto  threshold  would  be 
fine  for  verifying  the  concepts  and  algorithms  presented  in 
this  paper.  Table  4  lists  the  various  particles  found  in 
calculating  the  areas  and  centroids  for  each  image 
examined.  Particle  analysis  does  not  require  masking  so  no 
masks  are  indicated.  The  label  column  describes  the  image 
file  which  also  presents  the  source  positioning  information 
for  that  image.  For  example  the  label 

"upOdiam_ldiam_pos_x..."  indicates  that  the  source  light  was 
positioned  1  diameter  off  of  center  in  the  positive  X 
direction  and  was  not  moved  vertically,  etc.  Figures  74 
and  75  present  the  corresponding  coordinate  axes  for  both 
the  source  and  the  thruster  as  used  in  Table  4 . 


Label 

Area 

X 

Y 

Ma  j  or 

Minor 

Angle 

1 

on 

focus 

pt5secs 

4  aug  2005 

4343 

276.746 

29.401 

182 . 604 

30.282 

166.102 

2 

on 

focus 

pt5secs 

4  aug  2005 

16891 

292 . 18 

139.54 

160.747 

133.908 

134.72 

3 

on 

focus 

pt5secs 

4  aug  2005 

1175 

304.546 

429.943 

109.479 

13.665 

14.483 

97 


4 

upOdiam 

ldiam  pos  x  pt5secs  4  aug  2005 

6187 

185.416 

35.906 

183.007 

43.066 

11 . 87 

5 

upOdiam 

ldiam  pos  x  pt5secs  4  aug  2005 

27247 

258.621 

161 . 973 

218.615 

158 . 964 

120.044 

6 

upOdiam 

ldiam  pos  x  pt5secs  4  aug  2005 

2705 

203.531 

424 . 335 

166.181 

20.733 

164 . 91 

7 

upOdiam 

ldiam  neg  x  pt5secs  4  aug  2005 

3978 

337.551 

55.936 

173.186 

29.253 

147.823 

8 

upOdiam 

lpt5inch  neg  x  pt5secs  4  aug  2005 

5127 

314.891 

44.521 

189.146 

34.519 

154.33 

9 

upOdiam 

lpt5inch  neg  x  pt5secs  4  aug  2005 

6695 

334.321 

130.899 

117.068 

72.848 

133.903 

10 

upOdiam 

lpt5inch  pos  x  pt5secs  4  aug  2005 

5195 

178.505 

38.044 

168.821 

39.188 

13.164 

11 

upOdiam 

lpt5inch  pos  x  pt5secs  4  aug  2005 

28812 

253.873 

178 . 983 

225.809 

162 . 904 

116.247 

12 

upOdiam 

lpt5inch  pos  x  pt5secs  4  aug  2005 

3686 

187.472 

418 . 152 

188.46 

24 . 916 

161 . 517 

13 

upOdiam 

2diam  neg  x  pt5secs  4  aug  2005 

4431 

375.153 

86.678 

182.028 

30.994 

135.517 

14 

upOdiam 

3diam  neg  x  pt5secs  4  aug  2005 

4091 

383.473 

93.188 

175.383 

29.7 

133.809 

15 

upOdiam 

2diam  pos  x  pt5secs  4  aug  2005 

280 

182.303 

10.264 

30.403 

11.894 

13.253 

16 

upOdiam 

2diam  pos  x  pt5secs  4  aug  2005 

3431 

69.194 

94 . 311 

150.213 

29.099 

52.783 

17 

upOdiam 

2diam  pos  x  pt5secs  4  aug  2005 

22244 

203.33 

195.177 

224.069 

127 . 143 

121.013 

18 

upOdiam 

2diam  pos  x  pt5secs  4  aug  2005 

5181 

86.968 

354 . 649 

210 . 627 

31.325 

138.272 

19 

upOdiam 

3diam  pos  x  pt5secs  4  aug  2005 

1447 

45.048 

126.181 

78 . 912 

23.347 

65.158 

20 

upOdiam 

3diam  pos  x  pt5secs  4  aug  2005 

2634 

28.036 

295.441 

139.366 

24 .064 

119.143 

21 

upOdiam 

3diam  pos  x  pt5secs  4  aug  2005 

211 

176.264 

252.713 

24.213 

11 . 358 

113.949 

Table  4  Images  and  Area  Calculations  with  Centroids  and 
Ellipse  parameters . 


The  columns  in  Table  4  for  major,  minor  and  angle  are 
the  parameters  for  ellipses  that  approximate  the  areas  of 
light  in  the  images.  These  parameters  were  taken  for 
future  analyses  that  involve  using  the  major  axis  and  minor 
axis  values  to  determine  shrinkage  or  elongation  of  the 
light  areas  that  may  be  useful  to  increase  the 
effectiveness  of  the  area  moment  method  for  determining 
focal  spot  location.  The  angle  parameter  is  related  to  the 
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angle  between  the  major  axis  of  each  ellipse  and  a  line 
parallel  to  the  horizontal  or  X  axis  in  the  image.  The 
angle  and  the  length  measurements  could  be  combined  into 
the  method  of  area  moments  to  enhance  the  focal  spot 
location  determination. 


Figure  80  On-Focus  Image  For  Area  Measurements 

Figure  80  presents  the  situation  in  the  process  after 
the  image  has  been  cropped  and  a  threshold  has  been 
applied.  The  next  step  in  the  process  is  shown  in  Figure 
81  when  ImageJ  has  calculated  the  individual  areas  based  on 
a  minimum  size  of  100  pixels.  The  100  pixel  limit  came 
from  an  approximate  value  of  2000  pixels  that  would  be  in  a 
"tube  blob"  area  calculation  for  a  fully  lit  concentrator 
on-focus.  The  100  pixel  area  represented  3-5%  of  full 
light,  so  that  we  looked  for  light  areas  representing  95% 
fully  illuminated  concentrator.  The  values  for  the 
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centroids  and  areas  for  the  image  shown  in  Figure  81  are 
shown  in  table  1  with  label 

on_focus_pt5secs_4_aug_2005 .  Calculations  for  the  image: 
up0diam_ldiam_pos_x_pt5secs_4_aug_2005  gave  an  average  X 
value  of  215.856  and  average  Y  value  of  207.405.  The 
method  then  uses  either  the  centroid  of  the  central  "plate 
blob"  of  light  or  the  centroid  of  the  on-focus  image  to 
generate  the  difference  signal  used  for  concentrator 
control.  For  the  "tube  blob"  analysis,  the  difference 
calculation  gives  an  X  difference  of  -73.301  and  a  Y 
difference  of  -7.777.  By  subtracting  the  focal  spot 
centroids  from  the  centroid  of  the  "blob"  a  direction  of 
movement  towards  focus  is  generated.  If  the  on-focus  X  and 
Y  centroids  are  retained  from  that  image  a  larger 
displacement  value  may  be  generated  for  the  control  values 
that  move  the  focal  spot  towards  focus.  When  there  is  no 
central  "plate  blob"  of  light  in  the  image,  the  centroids 
for  X  and  Y  from  the  on-focus  image  have  to  be  used  for 
delta  calculations.  These  delta  calculations  would  then  be 
used  to  generate  control  values  for  the  hexapod  controller. 
Looking  at  Figure  82,  it  can  be  seen  that  this  method  works 
down  to  at  least  1.5  inches  off  of  focus  for  the  solar 
simulator.  That  number  is  associated  with  0.15  degrees  of 
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offset  for  the  test  setup,  which  is  very  close  to  the 


requirement  that  we  control  to  0.1  degree. 


Figure  81  Resulting  Particles  From  80,  Blobs  1  and  3  are  Concentrator 
Images  on  Tubes 


Centroid  Results 


Label 

1  on_focus_pt5secs_4_aug_2005 

2  on_focus_p  t5secs_4_aug_2005 

3  on_focus  j  t5secs_4_aug_2005 

Averages 


Area  X  Y 

4343  276.746  29401 

16891  292.18  13954 

1175  304546  429543 


Major  Minor  Angie 

182501  30282  166.102 

160.747  133508  134.72 

109479  13565  14483 


291.157  199.628 


4  up0diam_ldiam_pos_x_pt5secs_4_aug_2005 

5  up0diam_ldiam_pos_x_pt5secs_4_aug_2005 

6  upO  diam_l  diam_pos_x_pt5secs_4_aug_20Q5 


6187  185416  35506  183507  43566 

27247  258521  161573  218515  158564 

2705  203531  424335  166.181  20.733 


1137 

120544 

16451 


Averages 

Differences  (On  -  Idiam) 


215.856  207405 

73.301  -7.777 


Figure  82  Result  of  Centroid  Calculation 


Summary  Experiment  4: 
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Both  methods  work  for  locating  the  focal  spot  on  the 
sensor.  Correlation  using  only  Equation  15  fails  to  work 
below  2-3  diameters  of  misalignment.  Test  images  using 
high  pass  filtering  to  approximate  phase  only  correlation 
was  begun.  This  method  was  not  fully  implemented. 

Area  moments  did  better  below  2-3  diameters  of 
misalignment.  Figure  82  shows  the  complete  calculation  for 
one  set  of  images  used  in  this  study.  Good  focal  spot 
tracking  was  seen  from  5  diameters  of  misalignment  down  to 
1.5  inch  of  misalignment.  The  1.5  inch  gives  about  0.15 
degree  of  misalignment,  which  is  very  close  to  our  original 
requirement . 

Summary  and  Conclusions  of  Dissertation: 

A  system  for  locating  and  controlling  the  focal  spot 
of  a  solar  concentrator  has  been  developed  that  capitalizes 
on  in-situ  thruster  components  and  non-invasive 
measurements.  Compared  to  the  state  of  the  art  sensors 
described  in  section  3,  the  new  system  provides  enhanced 
capability  without  the  need  for  a  cavity  or  some  other 
intrusive  method  of  sensing  the  sunlight.  Also,  the  new 
system  is  not  limited  by  temperature  considerations  as  many 
of  the  sensors  in  section  3.  The  new  system  consists  of  a 
conical  shaped  absorber  doing  double  duty  as  the  heat 
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transfer  device  and  as  the  primary  sensor  that  is  imaged  by 
a  CCD  camera  as  the  measurement  device.  Both  of  these 
considerations  show  that  the  sensor  developed  in  this 
dissertation  is  superior  to  state  of  the  art  sensors  since 
none  of  the  sensors  in  section  3  can  provide  focal  spot 
information  for  continuous  tracking  of  the  Sun.  Wave  front 
sensor  techniques  are  introduced  and  modified  to  apply  to 
the  solar  concentrator  control  problem. 

Two  methods  of  tracking  the  focal  spot  were  then 
developed  by  modifying  the  wave  front  sensor  method.  One 
of  the  methods  developed  used  a  correlation  process  and 
masking  to  locate  the  focal  spot.  The  correlation  process 
was  shown  to  work  down  to  2-3  Sun  diameters  of  Sun 
misalignment.  The  other  method  called  area  moments  is 
developed  by  modifying  the  lenslet  spots  distance 
measurement  to  looking  at  areas  of  light  showing  in  the 
image.  By  calculating  an  average  area  centroid,  the 
location  of  the  focal  spot  can  be  determined  and  in  turn 
the  direction  to  move  the  focal  spot  can  be  determined. 

The  area  moment  method  was  shown  to  be  useful  down  to  an 
angle  of  misalignment  of  0.15  degrees.  This  angle  of  0.15 
degrees  is  very  close  to  the  arbitrary  requirement  of  0.1 
degrees . 
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An  area  for  future  work  is  possibly  a  combination  of 
both  methods  to  overcome  the  limitations  of  the  correlation 
method  alone.  Finally,  data  should  be  taken  below  1.5 
inches  of  misalignment  to  verify  that  the  method  can  get  to 
the  . 1  degree  of  misalignment  and  improve  the  control  down 
to  0.02  degrees  of  misalignment  where  the  thruster 
acceleration  could  be  utilized  in  the  optimization  of 
concentrator  location  using  a  combination  of  correlation 
and  area  moments  or  a  refined  area  moment  method. 


Appendix  A:  Section  4.2  Images 

Various  Photos  from  the  experimental  setup 


Figure  83  A1  CCD  Camera  Used  ST-6  From  SBIG 
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Figure  84  A2  Camera  Location  On  Strut  For  Concentrator 


Figure  85  A3  Concentrator  In  Torus  With  Camera  "on  Sun" 
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Figure  86  A4  Light  Source  With  Shielding 


Figure  87  A5  Light  Source  And  Scissor  Jack 
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Figure  88  A6  Thruster  Holding  Table 


Figure  89  A7  Absorber  and  Secondary  Concentrator 
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Appendix  B:  Java  Program  ConcentratorV3 

import  java.io.*; 
import  java. math.*; 
import  j avax . media . j 3d .* ; 
import  j avax . vecmath . * ; 

public  class  concentratorv3 


public  static  void  main (String [ ]  args)  throws  IOException 


System . out.println ("Hello  there ! " )  ; 

//  file  io  information  here 
InputStreamReader  reader  =  new 
InputStreamReader (System. in) ; 

Buf f eredReader  input  =  new  BufferedReader (reader ) ; 

System . out . print ( "Enter  a  filename  to  store  results:  "); 
String  fn  =  input.readLineO; 

System . out . println ( "This  is  the  file  name  you  gave  me:  "+ 

fn)  ; 


//System . out . println ( "Enter  a  filename  for  eij  output:  "); 
//String  fn2  =  input.readLineO; 

//System . out . println ( "This  is  the  second  file  name  you  gave 
me :  "  +  f n2 ) ; 

//ends  for  input  of  filename. 

int  N  =  20;  //number  of  points  in  grid  =  (2*N-1) * (2*N-1) 
double  R  =  7.0;  //radius  of  projected  circle 
//double  f  =  7.28;  //  focal  point  location 

double  f  =  7.28; 

double  theta  cone  =  30;  //half  angle  of  light  cone (degrees ) 
double  phi_cone  =90;  //  tilt  angle  of  light 

cone . (degrees ) 

int  li  =  0; 

double  mzh  =  1.52e9;  //  sun  distance  change  to  correct 

val . 


Vector4d[]  intensity  =  new  Vector4d [ (2 *N-1 ) * (2 *N-1 ) ] ; 
//  image  of  light  in  concentrator 
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double  x  =0.0;  double  y  =  0.0;  double  thz  =  0.48;  double 

phz  =  170.0; 

//  x,  y  are  the  deltas  in  phi  and  theta_c.  thz  theta_zhi 
and  phz  phi_zhi 

boolean!]  xcirc  =  new  boolean [ (2*N-1 )* (2*N-1 )]; //  vector  to 
store  logicals  defining  xs  in  circle 

Vector3d[]  xl  =  new  Vector3d [ (2*N-1 ) * (2 *N-1 ) ] ; //  xl  vector 
Vector3d[]  xtemp  =  new  Vector3d [ (2*N-1 ) * (2*N-1 ) ] ; //  temp 

vector 

Vector3d  s  =  new  Vector3d (0 . 0, 0 . 0, f ) ;  //focus  vector 

Vector3d  st  =  new  Vector3d (0 . 0, 0 . 0, f ) ;  //  copy  of  focus 

vector 

Vector3d  sta  =  new  Vector3d (0 . 0, 0 . 0, 0 . 01) ; //vector  to  move 


eye 

vectors 


Vector3d[]  eij  =  new  Vector3d [ (2*N-1 ) * (2*N-1 ) ] ; //  pointing 


Vector3d  r,T,zhi;  //  mirror  to  Sun  vector,  reflected 
vector,  origin  to  sun  vector 

Vector3d[]  TT  =  new  Vector3d [ (2*N-1 ) * (2 *N-1 ) ] ; 

Vector3d[]  image  =  new  Vector3d [ (2*N-1 ) * (2*N-1 ) ] ; //  image 
of  sun  on  focal  plane 

Vector3d[]  eijtemp  =  new  Vector3d [ (2*N-1 ) * (2 *N- 
1 )]; //variable  to  hold  original  eij  vectors, 
double  dth  = 

(2 . 0*theta_cone* (Math. PI/180.0) ) / (2.0* (double ) N-l ) ; 
double  dphi 

= (2 . 0*theta_cone* (Math. PI/180.0) ) / (2.0* (double ) N-l ) ; 
double  zx  = 

mzh*Math .sin ( thz*Math .PI/180.0) *Math .cos (phz*Math .PI/180.0) ;// vector 
components  for  zhi 

double  zy  = 

mzh*Math .sin ( thz*Math . PI / 1 80 . 0 ) *Math .sin (phz*Math . PI/ 1 80 . 0 )  ; 

double  zz  =  mzh*Math . cos (thz*Math . PI/180 . 0 ) ; 
double  costhet  =  Math . cos ( theta_cone*Math . PI / 180 . 0 );/ /used 
to  make  sure  eij  vectors  are  all  inside  the  cone  of  light 

zhi  =  new  Vector3d ( zx, zy, z z ) ;  //vector  from  origin  to  sun 
System . out . println (" ZHI  =  "  +  zhi); 

Vector3d  eoo  =  new 

Vector 3d ( 0 , Math .sin (phi_cone*Math . PI/ 180 . 0 ) , Math .cos (phi_cone*Math . PI/ 1 
80.0));//  e  zero  zero  vector 

System . out . println ( "EOO  =  "  +  eoo)  ; 

Vector3d[]  incid  =  new  Vector3d [ (2 *N-1 ) * (2*N-1 ) ] ; // 
incident  vectors  to  mirrors 

Vector3d[]  norm  =  new  Vector3d [ (2*N-1 )* (2*N-1 )]; //normal 


vectors 

Vector3d[]  nn  =  new  Vector3d [ (2*N-1 ) * (2 *N-1 ) ] ; //  copy  of 
norm  values  to  store 

double  cos_sun_t  =  Math . cos ( 0 . 5*Math . PI/ 1 80 . 0 ); //cosine  of 
angle  of  sun  disk 

for(int  i  =  0;  i  <  2*N-l;i++) 

{ 

for(int  j =0 ; j <2*N-1 ; j ++) 

{ 


(0.0*Math. PI/180) ; 


x  =  ( (double) i  -  (double) N+  1.0)*dth  + 

y  =  ( (double) j  -  (double) N  +1.0)*dphi; 
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//  don't  forget  to  change  the  x,  y  below  if  you 

change  them  here 

ei j [ li ]  =  new 

Vector 3d (Math . tan (x) , 0 , Math . tan (y) ) ; 

eij  [li]  .add(eoo) ; 
ei j [ li ] . normalize ( ) ; 

//System . out . println ( "eij [ "  +  li  +  "]  =  "  + 

eij  [li]); 

//System . out . println ( "x  =  "  +  x  +  "y  =  "  +  y  + 

"  and  zhi  =  "  +  zhi  +  "li  =  "  +  li) ; 

li++;  //  increment  matrix  index 


//  Start  of  intercept  and  mirror  element  calculation 

li  =  0; 

st . add (sta) ; 

st . scale  (-1.0) ; 

double  d  =  0; 

for(int  i  =  0 ; i<2*N-l ; i++) 

{ 

for(int  j  =  0 ; j <2*N-1 ; j ++) 

{ 

double  a  = 

(eij  [li]  . x) * (eij [li]  . x)  +  (eij [li]  . y) * (eij [li]  .y) ; 

double  b  =  -4 . 0*f*eij [li] . z; 

double  c  =  -4.0*f*f;  //  a,  b,  c  from  ax^2  + 

bx  +  c  =  0 

d  =  b*b-4.0*a*c;  //discriminant 
double  tl,t2,t=0; 

if (a  ==  0.0)  System . out . println ( "x  =  "  +  -c/b) ; 
if (d  ==  0.0)  System . out . println ( "x  =  "  +  - 

b/  (2 . 0*a) ) ; 

if (d  >  0.0) 

{  //System . out . println ( "Two  real  roots") ; 

tl  =  (-b  +  Math . sqrt (d) ) / (2 . 0*a) ; 
t2  =  (-b  -  Math  .  sqrt  (d)  )  /  (2 . 0*a)  ; 
if(tl  >  0.0  &&  t2  <  0.0)  t  =  tl; 
if(t2  >  0.0  &&  tl  <  0.0)  t  =  t2; 

} 

if (d  <  0.0) { 

System .out. println (" Imaginary  roots" ) ; 
System . exit  ( 0 ) ; 

} 

//  now  find  x's  which  are  the  mirror  element 

position  vectors 

xl  [  li ]  =  new 

Vector3d ( t*ei j [ li ] . x, t*ei j [ li ] . y, t*ei j [ li ] . z  +  f ) ; 

//  now  find  normal  vector  to  each  xl 
nn[li]  =  new  Vector3d ( 0 , 0 , 0 ) ; 
norm[li]  =  new  Vector3d(- 
xl  [  li]  . x/ (2 . 0*f ) / (Math . sqrt (xl [li]  . z/f  +  1.0)),- 
xl  [ li ]  . y/ (2 . 0*f ) / (Math . sqrt (xl [ li ]. z/f  + 

1 . 0) ) , 1 .0/ (Math. sqrt (xl  [li]  .z/f  +1.0))); 

//already  normalized 

// 
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norm [ li ] ) ; 
xl [li] )  ; 


vectors 

calc  eij  vectors  only 


incid [ li ] ) ; 


nn [ li ]  . add (norm [ li ] )  ; 

//System . out . println ( "norm [ "  +  li  +  "]=  "  + 

//System . out . println ( "xl [ "  +  li  +  "]  =  "  + 

incid[li]  =  new  Vector3d ( 0 , 0 , 0 ) ; 

incid [ li ]  . add (xl [ li ] )  ; 

eijtemp[li]  =  new  Vector3d ( 0 , 0 , 0 )  ; 

eij temp [li] . add (eij [li] ); //  store  original  eij 

//changed  way  of  calculating  incid  as  original 

//eij [li] .scale (-1.0) ; 

//incid[li]  . add (eij [ li ]  )  ; 
incid[li] .add(st) ; 

//System . out . println (" incid [ "  +  li  +  "]  =  "  + 
li++;  //increment  matrix  index 


//  now  we  have  arrays  of  Xs,  norms,  eyes,  and  incident  rays 
//  using  the  law  of  reflection,  will  calculate  t,  reflected 

ray 

//  T=incid[li]-  2* (norm [li] . dot (incid [ li] )) *n 
//  if  T.R  /  | T | * | R |  >=  cos  theta_sun  then  "hot" 

//  if  T.R  /  | T | * | R |  <  cos  theta_sun  then  "cold" 

//  if  "hot"  store  power  of  sun  portion  into  intensity 
//  generate  array  of  decision  values  for  circular  shadow 
li  =0; 

for(int  i  =  0 ; i<2*N-l ; i++) 

{ 

for(int  j  =  0 ; j <2*N— 1 ; j ++) 

{ 

if (xl [li ] . x*xl [ li ] . x+ (xl [ li ] .  y- 
14 . 81) * (xl [li] .y-14 . 81)  <=  R*R) 

xcirc[li]  =  true; 

else 

xcirc[li]  =  false; 

1  i  +  + ; 


li  =  0; 

double  talpha  =  0.0; 
double  dottn  =  0; 

//open  file 
try 
{ 

FileWriter  out  =  new  FileWriter (fn) ; 

Buf f eredWriter  brOut  =  new  Buf feredWriter (out) ; 
for(int  i  =  0;  i  <  2*N-l;i++) 

{ 

for(int  j  =  0;  j<2*N-l;j++) 

{ 

x  =  ( (double) i  -  (double) N+  1.0)*dth  + 

(00. 0*Math. PI/180)  ; 
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y  =  ( (double) j  -  (double) N  +1.0)*dphi; 

//System . out . println ( "x  =  "  +  x  +  "  y  =  "  +  y)  ; 
//System . out . println (" li  =  "  +  li); 
r  =  new  Vector3d ( 0 , 0 , 0 ) ; 

T  =  new  Vector3d ( 0 , 0 , 0 ) ; 
intensity [ li ]  =  new  Vector4d ( 0 , 0 , 0 , 0 ) ; 
image [li]  =  new  Vector3d ( 0 , 0 , 0 ) ; 
xtemp [li]  =  new  Vector3d ( 0 , 0 , 0 ) ; 

TT[li]  =  new  Vector3d ( 0 , 0 , 0 ) ; 
xtemp [ li ] . add (xl [ li ] ) ; 
xtemp [li]  . scale  (-1.0); 
r . add ( zhi ) ; 

r . add (xtemp [li ]) ;  //calculate  ray  from  mirror 

r . normalize ( ) ; 

//incid[li]  .scale  (-1.0) ; 

T . add ( incid [ li ] ) ; 

//System . out . println ( "T  =  "  +  T  +  "at  li  =  "  + 

//System . out . println ( "eij  =  "  +  eij[li]  +  "at 

dottn  =  -2 . 0* (norm [ li] . dot (incid [ li] ) ) ; 
norm [ li ] . scale (dottn) ; 

T . add (norm [ li ] ) ; 

TT [li] . add (T) ; 

//System . out . println ( "T  =  "  +  T)  ; 

//System . out . println ( "R  =  "  +  r) ; 
talpha  =  T . dot (r) / (T . length () *r . length ()) ; 
//System . out . println  (  "  i  =  "  +  i  +  "j  =  "  +  j); 
//System . out . println ( "talpha  =  "  +  talpha); 
if (xcirc  [ li ] ) 

{ 

if ( (eij  temp [li]  . dot (eoo) / (eij  temp [ li ]  . length ( ) *eoo . length ( ) )  >= 
costhet)  &&  talpha  >=  cos_sun_t) 

{ 

intensity [ li ]  =  new 

Vector4d(xl[li]  . x, xl [ li ]  . y, xl [li ]  . z , 1 . 0 )  ; 

image [li]  =  new  Vector3d (x, y, 1 . 0 ) ; 

} 

else 

{  intensity [ li ]  =  new 

Vector4d(xl[li] . x, xl [ li ] . y, xl [li ] . z , 0 . 0 ) ; 

image [li]  =  new  Vector3d (x, y, 0 . 0 ) ; 

} 

} 

else 

{  intensity [ li ]  =  new 

Vector4d(xl[li]  . x, xl [ li ]  . y, xl [li ]  . z  ,  0 . 0 )  ; 

image [li]  =  new  Vector3d (x, y, 0 . 0 ) ; 

} 

//System . out . println (" intensity [ "  +  li  +  "]  =  " 

+  intensity[li]); 

//System . out . println (" image [ "  +  li  +  "]  =  "  + 

image [ li ] ) ; 

//write  file 


to  sun 

li)  ; 

li  =  "  +  li) ; 
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//brOut . write (intensity [ li] . x  +  "\t"  + 
intensity [ li ]. y  +  "\t"  +  intensity [ li ]. z  +  "\t"  +  intensity [ li ]. w  + 

"  \  n" )  ; 

brOut . write ( image [ li ]. x  +  "\t"  +  image [ 1 i ] . y  + 
"\t"  +  image [li].z+  "\t"  +  li  +  "\t"  +  "\n"); 

1  i  +  + ; 


//close  file 
brOut . close ( ) ; 
out . close  ( ) ; 

} 

catch  (IOException  e)  { 


} 

//File  write  to  save  incid,  T,  and  norm  values 
/*li  =  0; 


try 

{ 

FileWriter  out  =  new  FileWriter ( fn2 ) ; 

Buf f eredWriter  brOut  =  new  Buf feredWriter (out) ; 
for(int  i  =  0;  i  <  2*N-l;i++) 

{ 

for(int  j  =  0;  j<2*N-l;j++) 

{ 

brOut . write ( incid [ li ]. x  +  "\t"  +  incid[li].y  + 
"\t"  +  incid [ li ] . z+  "\t"  +  xl[li].x  +"\t"  +  xl[li].y  +  "\t"  +xl(li].z 
+"\t"  +  nn[li] ,x+"\t"+  nn[li] . y+  "\t"  +  nn[li] .z  +  "\t"  +li  +  "\t"  + 

"  \  n" )  ; 

1  i  +  + ; 


//close  file 
brOut . close ( ) ; 
out . close  ( ) ; 

} 

catch  (IOException  e)  { 
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